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Work  under  this  Contract  has  been  initiated  in  several 
directions.  In  one  major  and  one  minor  projects  there  has  been 
sufficient  progress  to  warrant  publication.  These  are  discussed 
below  and  detailed  accounts  are  enclosed.  Two  additional  projects 
which  appear  promising  require  additional  work  and  are  only  ment- 
ioned here  (under  (c)  below)  as  a preliminary  account. 

(a)  A theoretical  framework  for  the  interpretation  of  structural 
effects  in  reactive  molecular  collisions  has  been  formulated. 

The  approach  has  been  presented  in  a number  of  seminars  and  at  the 
NBS  meeting  'Current  Status  of  Kinetics  of  Elementary  Gas  Reactions  ' 
A detailed  account,  entitled  'Reactivity,  Energy  and  Structure', 
has  been  prepared  for  publication  (in  the  Journal  of  Physical 
Chemistry),  and  is  enclosed  as  part  of  this  report.  The  formalism 
: developed  therein  is  primarily  thermodynamic-like  in  character  and 

hence  fails  to  provide  a mechanistic  interpretation.  A first 
attempt  to  provide  one  (section  7)  was  restricted  to  reactions  on 
a given  electronic  surface.  Future  work  will  examine  a more  refined 

approach  aimed  at  an  extension  to  processes  including  electronic 

I 

excitation . 

^ (b)  An  exact  information  theoretic  solution  to  the  problem  of  a 
collinear  collision  of  an  atom  with  a Horse  type  (i.e.  anharmonic) 
oscillator  has  been  provided.  An  account  has  been  prepared  for 
publication.  The  Work  demonstrates  a novel  technique  for  the 
determination  of  constraints  in  problems  which  at  first  glance 
appear  insoluble. 

(c)  Preliminary  work  has  been  carried  out  on  the  question  of  the 
exact  definition  of  a 'prior'  expectation  for  the  overall  magnitude 

» \ " ’ 
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of  the  rate  constant.  Such  a definition  is  particularly  relevant 
to  electronic  excitation  processes  because  of  the  large  differences 
in  the  available  energy  that  can  occur  between  the  different 
reaction  products.  It  appears  that  a natural  definition  free  from 
ambiguity  is  indeid\feasible . The  basis  for  these  developments  is 

the  appendix  to  the  manuscript  'Reactivity,  Energy  and  Structure' 

\ 

\ 

which  is  enclosed.  A Related  development  has  been  the  search  for 
a formalism  that  will  explicitly  recognise  that  often  one  is  only 
able  to  make  some  general  statements  about  the  Hamiltonian  of  the 
system  but  cannot  specify  it  in  detail.  Preliminary  results  are 
encouraging  and  numerical  experiments  (in  collaboration  with  W.  S. 
Warren,  S.  Mukamel  and  A.  Pines)  are  in  progress. 


\t 
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Reactivity,  Energy  and  Structure 


R.  D.  Levine 
Department  of  Chemistry 
The  University  of  California 
Berkeley,  California  94720 


Abstract 

The  reaction  rate  constant  is  expressed  as  Zexp(-G _/RT).  Z 

a 

is  the  binary  collision  frequency.  G_ , the  free  energy  of  activation 

Cl 

is  shown  to  be  the  difference  between  the  free  energy  of  the  reactive 
reactants  and  the  free  energy  of  all  reactants.  The  results  are 
derived  from  both  a statistical  mechanical  and  a collision  theoretic 
point  of  view.  While  the  later  is  more  suitable  for  an  ab-initio 
computation  of  the  reaction  rate,  it  is  the  former  that  lends  itself 
to  the  search  of  systematics  and  of  correlations  and  to  compaction 
of  data.  Different  thermodynamic-like  routes  to  the  chatacterisat- 
ion  of  G are  thus  explored.  The  two  most  promising  ones  appear 
to  be  the  use  of  thermodynamic  type  cycles  and  the  changes  of 
dependent  variables  using  the  Legendre  transform  technique.  The 
dependence  of  G=  on  AG°,  the  standard  free  energy  change  in  the 

a 

reaction,  is  examined  from  the  later  point  of  view.  It  is  shown 
that  one  can  rigorously  express  this  dependence  as  G = aAG°  + 

d 

G°M(a).  Here  a is  the  Brpnsted  slope,  a = -31nk(T)/3(AG°/RT) , G° 

d d 

is  independent  of  AG°  and  M(a),  the  Legendre  transform  of  G , is 

d 

< 


a function  only  of  a.  For  small  changes  in  AG°,  the  general  result 
reduces  to  the  familiar  'linear'  free  energy  relation  6G=  = a6AG°. 

a 

It  is  concluded  from  general  considerations  that  M(o)  is  a symmetric, 
convex  function  of  a and  hence  that  a is  a monotonically  increasing 
function  of  AG°.  Experimental  data  appear  to  conform  well  to  the 
form  a = 1/Cl  + exp(-AG°/G°) ] . A simple  interpretation  of  the 

CL 

AG°  dependence  of  G , based  on  an  interpolation  of  the  free  energy 
from  that  of  the  reagents  to  that  of  the  products  is  offered. 


Introduction 

Thermodynamics  is  distinguished  not  only  by  the  systematic 

statement  of  its  principles  but  also  by  the  considerable  compactness 

and  order  in  its  data  base.  This  paper  explores  the  question 

whether  thermodynamic-like  methods  can  be  formulated  towards 

2 3 

achieving  similar  goals  in  chemical  kinetics.  ’ Such  methods  have 

already  been  explored  for  relative  rate  constants,  with  special 

4 

emphasis  on  energy  disposal  and  energy  requirements.  Here  however 

we  are  particularly  concerned  with  examining  the  magnitude  of  the 

overall  rate  constant  and  its  variation  both  with  energy  (e.g. 

temperature)  and  with  structural  modifications  of  the  reactants. 

5 — 8 

The  search  for  such  systematics  preceeded  even  transition  state 
3 9 2 3 10—2  8 

theory  * , and  has  continued  * * ever  since.  The  point  of 

view  adopted  here  will  however  be  somewhat  different  and  will 

5 6 

follow  a line  of  reasoning  started  by  Tolman  and  La  Mer  . To 

reduce  the  formal  results  to  concrete  terms,  the  operational 

7 

definition  of  the  Br^nsted  exponent  and  the  idea  of  chemical 
inertia  due  to  Evans  and  Polanyi^  are  employed.  In  a differ- 
ent vein,  it  is  explored  whether  the  familiar  concept  of  a 
thermodynamic  cycle  can  be  adapted  for  applications  in  chemical 
kinetics . 

Specifically,  this  paper  offers  a derivation  and  simple 
applications  of  the  concept  of  the  free  energy  of  activation. 

The  derivation  is  based  on  the  introduction  of  the  'distribution 
of  reactive  reactants'^  using  either  a collision-theoretic  or 
a statistical  mechanics  point  of  view.  The  free  energy  of 


activation  is  then  shown  to  be  the  free  energy  of  the  distribution 
of  reactive  reactants.  Applications  include  both  rigorous  results 


(e.g.  cycles)  and  approximations  based  on  the  idea  of  a Br^nsted 

slope.  The  type  of  kinetic  data  that  can  be  used  to  advantage  to 

test  the  theory  and  to  guide  future  work  is  noted. 

Reactivity,  or  the  fraction  of  collisions  that  lead  to  reaction, 

is  discussed  from  a thermodynamic-like  point  of  view  in  section  II. 

4 29-31 

It  was  previously  argued  ’ that  one  can  introduce  thermodynamic 

functions  and  hence  cycles  for  systems  not  in  equilibrium.  Section 

III  explores  the  possibility  of  comparing  reactivities  in  this 

fashion.  Other  theoretical  routes  to  the  determination  of  the 

reactivity  are  examined  in  section  IV.  Finally,  section  V provides 

a discussion  of  the  reactivity  from  a collision-theoretic  point  of 

view,  leading  to  the  expression 

k (T)  = Zp  (T)  (1.1) 

r r 

for  the  reaction  rate  where  Z is  the  collision  frequency  and  Pr(T) 
is  the  reactivity.  An  appendix  relates  this  result  to  transition 
state  theory,  by  showing  that 

Jr(T)  = kTST(T)K(T).  (1.2) 

Here  kTgT(T)  is  the  transition  state  theory  rate  constant  and  the 

transmission  coefficient  k(T),  k(T)<  1,  is  the  reactivity  of  the 

’transition  state'  i.e.  the  fraction  of  collisions  that  reach  the 


transition  state  and  proceed  to  form  products  Turned  around,  the 
expression  kr(T)  = Zpr(T)  can  be  regarded  as  employing  transition 
state  theory  with  a very  very  early  transition  state  (so  that 
k_  (T)  = Z)  but  with  an  exact  account  of  the  transmission  coeffic- 

lu  1 
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Structure-reactivity  relations  are  discussed  in  sections  VI 
and  VII.  Starting  with  the  definition  of  the  Brynsted  slope  a 
a=-3lnkr(T)/3(AG°/RT)  as  the  change  in  the  rate  constant  with  the 
free  energy  change  of  the  reaction  it  is  shown  that  for  any 
Br0nsted  series  one  can  rigorously  express  kr(T)  as 

kr  = Zexp(-Ga/RT)  (1.3) 

Here 

G = otAG°  + G °M(a)  (1.4) 

where  Z is  independent  of  AG°.  M(a)  is  a function  of  a defined  by 
3M(ct)/3a=-AG°/G° 

3 

so  that  5G  = ct6AG°  and  G°  is  a parameter,  with  the  dimensions  of 

3 3 

energy,  which  is  characteristic  of  the  reaction  series.  The  mono- 
tonic increase  of  a,  and  hence  of  G , with  AG°  (i.e.  the  reactivity- 
selectivity  principle)  are  shown  to  be  valid  for  any  Brpnsted 
series.  Section  VII  examines  a simple  model  for  the  variation  of 
free  energy  enroute  from  reagents  to  products.  The  concept  of  the 
inertia  of  chemical  reactions,  due  to  Evans  and  Polanyi10  is 
shown  to  be  sufficient  to  imply  that  the  Brdnsted  slope  a can  be 

interpreted  as  the  location  of  the  free  energy  barrier  to  reaction 

12 

(the  Leffler  postulate  ) and  that  a shifts  towards  the  products 

o 13  32 

side  as  AG  increases  (the  Hammond  postulate  ;.  The  correlation 


of  barrier  location  with  energy  requirements  and  energy  disposal 
of  chemical  reactions  and  with  branching  ratios  is  also  discussed. 


II.  A Thermodynamic  Like  Approach 


The  purpose  of  this  section  is  to  offer  a derivation  and 
interpretation  of  the  expression 

p = exp (-G  /RT)  (2.1) 

r oL 

for  the  reactivity  in  terms  of  a free  energy  G . Since  (2.1)  will 
be  derived  as  an  identity,  valid  irrespective  of  the  detailed 
mechanism  of  the  collision,  it  is  of  merit  only  if  one  can  forge 
an  independent  route  to  G= . Attempts  in  this  direction  are  dis- 
cussed  in  the  subsequent  sections  of  this  paper  and  additional 
routes  are  under  active  examination. 

The  definition  of  reactivity  to  be  introduced  in  this  section 
is  based  on  the  following  physical  picture.  For  a bimolecular 


reaction,  consider  a large  number  of  independent  binary  collisions 

(as  in  the  case  in  a molecular  beam  or  a chemiluminescence  experi- 
33 

ments  ) . The  initial  states  of  the  colliding  molecules  are 

randomly  selected  from  a system  in  thermal  equilibrium  at  the 

temperature  T.  Not  all  such  collisions  necessarily  lead  to 

formation  of  products,  as  some  may  be  non  reactive.  Reactivity  is 

the  fraction  of  such  collisions  that  do  lead  to  reaction.  In  this 

section,  and  then  again  in  section  III  we  are  concerned  with  a 

thermodynamic- like  characterisation,  while  a corresponding 

collision- theoretic  analysis  is  provided  in  V.  To  make  sure  that 

the  discussion  retains  an  intuitive  flavor,  we  shall  use  concepts 

34  35 

appropriate  to  a quasiclassical  description  of  the  collision.  ' 


In  a computational  study,  the  independent  binary  collisions 
are  simulated  by  running  N classical  trajectories  (where  N is 


! 

5- 

a large  number) . The  trajectories  differ  in  their  initial  condit- 
ions and  these  are  chosen  to  represent  the  distribution  of  the 
34, 36 

reactants. 

Thermal  Reactants . To  simulate  a thermal  distribution  of  the  reagents 
one  needs  to  run  trajectories  at  different  total  energies.  Since 
only  a finite  number  of  trajectories  will  be  examined  one  needs  to 
represent  the  distribution  of  total  energy  as  a discrete  histogram, 
such  that  there  are  trajectories  with  the  (total)  energy  E^.  The 
average  energy  (per  trajectory)  is 

<E>=ENiEi/N.  (2.2) 

For  a system  at  thermal  equilibrium,  <E>  suffices  to  characterise 
the  distribution  of  energy  in  the  reactants.  We  show  this  explicitly, 
as  a preliminary  to  the  next  stage  of  the  argument. 

Let  there  be  gi  distinct  initial  state  of  the  energy  E^ . The 

4 

entropy  of  the  initial  ensemble  is  (in  units  of  the  gas  constant) . 

S =-Ep.lnp.  + Ep . S . (2.3) 

i i l i 1 1 

■ -Epiln(pi/gi) 

where  p^  = N^/N  and  = lng^.  If  p^  = g^/N  then  S = InN.  However, 
because  of  the  constraint  (2.2)  it  is  not  feasible  to  have  p^=g^/N. 

Rather,  one  needs  to  determine  the  p. 's  subject  to  <E>=EE.p.  and 

i 4 

to  E^p^  = 1 by  maximising  the  entropy.  The  result  is  well  known 

p±  = giexp(-Ei/RT)/Q  (2.4) 

where  Q insures  that  the  p^s  are  normalised, 

Q = Egiexp(-Ei/RT)  (2.5) 

and  the  magnitude  of  T is  adjusted  to  reproduce  the  given  magnitude 

of  E.  In  other  words,  T is  that  number  that  satisfies 

<E>=EE.g.exp(-E./RT) /Eg . exp (-E . /RT) =RT231nQ/3T.  (2.6) 

i 1 1 1 i 1 1 

J 

■ . A 


Reactive  Reactants.  The  equations  of  motion  allow  one  to  determine 
whether  a given  trajectory  is  reactive  or  not.  Hence,  if  Nr^  is  the 
number  of  reactive  trajectories  with  the  total  energy  E.  and  N . = 


Nr^  is  the  number  of  non  reactive  ones 


The  central  point  is  that  not  only  the  energy  but  also  other 


To  show  this,  it  is 


ies  with  the  energy  E 


"ri'  i'  ri'  "i*  rr  ri,ri 

has  the  interpretation  of  the  fraction  of  trajectories  from  the 
state  i which  are  reactive  and  will  be  denoted  g ./g. , so  that 


9i  " gri  + 9nri*  (J-13> 

It  is  important  to  note  that  even  if  we  refine  the  interpretation 

of  the  index  i,  up  to  and  including  i being  a complete  specification 

of  an  initial  quantum  state,  it  would  still  be  the  case  that 

0 < gr^  < 1 (rather  than  gri  being  either  unity  of  zero) . The 
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reason  is  that  a single  quantum  state  is  simulated  by  a group  of 
classical  trajectories  (which  differ  in  their  initial  classical 
conditions  and  while  every  single  trajectory  is  either  reactive 
or  not,  the  different  members  of  the  group  may  behave  differently) 


Reactivity.  The  resolution  (2.13)  implies  a similar  resolution 


of  the  partition  function 


0 = 0 + Q 

v vr  nr 


where 


Qr  « Egri  exp  (-E^/RT) 


(2.14) 


(2.15) 


and  similarly  for  Q . With  this  definition  of  Q„  and  the  results 
■*  nr  r 

prxri  = Pi9ri/'3i  and  (2.12)  we  immediately  obtain 


= » r*ri 

■ ?Pi(9ri/9i) 

1 

= Qr/Q 


(2.16) 


This  short  derivation  tends  to  stress  that  the  result  pr  = Qr/Q  is 
indeed  an  identity  but  avoids  a proper  interpretation  of  Qf.  A 
longer,  but  more  enlightening  route  is  as  follows:  in  terras  of  the 
reactivity,  the  expression  (2.3)  for  the  entropy  can  be  brought 


to  the  form 


S - -pr1npr-pnr1npnr+prsr+pnrs, 


where 


S = -lx  . lnx  .+Ex  .S 
r ^ ri  n i ri  n 


(2.17) 


(2.18) 


= -Zx  . In (x  ./g  , ) 
i ri  rx'  *ri 


with  Sri  = lngri . 
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Using  (2.17)  and  (2.18)  one  can  now  repeat  the  process  of 
searching  for  the  distribution  of  maximal  entropy  subject  to  the 
condition  of  given  mean  energy 


<E>  = p Ex  .E.  + p Ex  .E. 

*r^  ri  i nri  nri  *• 

and  normalisation 


(2.19) 


p + p =1 

*r  *nr 


(2.20) 


Exri  1/  2xnri  - 1. 


The  result  is 


(2.21) 


xri  = ^riexP (~Ei/RT) /Qr 
Qr  = Egriexp(-Ei/RT) 

pr  = exp(Sr-<E>r/RT)/Q 

Q = exp (S  ~<E>  /RT)  + exp (S  -<E>  /RT) 
r r nr  nr 

Substituting  (2.22)  in  (2.18)  we  have,  using  (2.15),  that 


(2.22) 

(2.15) 

(2.23) 

(2.24) 


sr  = lnQr  + <E>r/RT 

and  hence  derive  from  (2.23)  the  desired  results 

Pr  = Qr/Q 

Q = Q + Q 

r nr 


(2.25) 


(2.16) 

(2.14) 


For  a system  with  a given  value  of  the  mean  energy,  it  is 
the  partition  function  that  determines  pr.  In  terms  of  the  free 


energy 


Here 


Gr  = -RTlnQr 
pr  = exp [- (Gr~G) /RT] 
= exp [-G  /RT] 

a 

= exp [S  - E /RT] 

cl  a 


Ga  = Gr-G 


(2.26) 

(2.27) 


(2.28) 
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Ea  = <E>r  " <E>  (2.29) 

and 

Sa  = Sr  - S.  (2.30) 

Pending  further  discussion  below  it  should  already  be  stated  here 
that  Sr  is  not  the  information-theory  entropy  of  the  distribution 
of  the  reactive  reactants  (but  is  related  to  it  via  a Legendre 
transform) . 

Energy^  Entropy  and  Work.  To  complete  the  discussion  we  note  two 
additional  implications.  The  first  is  the  validity  of  a Gibbs- 
Helmholtz  like  theorem.  To  see  this  note  that  from  (2.10)  and  (2.26) 
<E>r  = RT231nQr/3T  (2.31) 

= -T23 (Gr/T)/3T 
and  hence  that 

Sr  = _3Gr/3T-  (2.32) 

In  a more  technical  terminology  we  can  write  these  identities  as 

Sr  = lnQr  + T23(lnQr)/3T  (2.33) 

. . 37 

showing  that  Sr  is  the  Legendre  transform  of  lnQr-  In  practice 

these  identities  imply  that  pr  has  a well-defined  temperature 

dependence 

-31np  /3 (1/T) =RT231np  /3T  = Ea  = <E>  -<E>  (2.34) 

It  is  important  to  note  (as  was  first  pointed  out  by  Tolman'*'2®) 
that  (2.34)  is  an  exact  relation  and  is  by  no  means  conditioned  on 
Efl  being  (nearly)  temperature  independent.  Whatever  the  temperature 
dependence  of  Ea  may  be,  the  relation  (2.34)  holds.  The  reason 
is  that  S„  satisfies  the  Gibbs-Helmholtz  like  reaction 
Sa  = -3G/3T 

3 cl 


(2.35) 


or,  equivalently. 


3S  /3T  = T_13e  /3T . (2.36) 

d cL 

which  is  invoked  in  proving  (2.34).  If  one  defines  C_  = 3E  /3T  as 

cl  cl 

the  difference  in  the  specific  heats  between  the  reactive  reactants 
and  the  reactants  then 

3S  /3T  - C /T  (2.37) 

cl  cl 

is  the  relation  first  noted  by  LaMer.6 

The  second  point  is  that  -G,  = G-G„  does  admit  of  an  inter- 

pretation  as  the  maximal  work  when  a mole  of  reactive  reactants 

(at  the  temperature  T)  is  brought  to  thermal  equilibrium  (at  the 

same  temperature) . The  general  approach  to  the  computation  of 

thermodynamic  functions  for  systems  in  disequilibrium  has  been 

29  31 

previously  discussed.  ' Here  we  only  need  the  final  result. 

If  x^  is  the  actual  fraction  of  systems  in  the  state  i,  and  pr^ 
is  the  thermal  equilibrium  fraction  of  systems  in  the  state  i at 
the  temperature  T,  then  the  maximal  work  available  in  an  isothermal 
process  at  the  temperature  T is,  on  statistical  mechanical  grounds30, 

RT?Xriln(xri/pri}  = "Ga  (2.38) 

where  the  equality  (2.1)  follows  when  one  notes  that  pri  = Prxr^* 

If  one  prefers  purely  thermodynamic  reasoning,  the  maximal  work 

39 

under  isothermal  conditions  is 

( <E>-<E>  ) -RT (S-S  ) = - (E  -RTS  ) ■ -G=  (2.39) 

I.  i ad  a 

Of  course,  the  simplicity  of  (2.39)  is  misleading  for  without  inform- 
ation theory  Sr  is  not  easily  interpreted. 
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Summary . From  a purely  algebraic  point  of  view  what  was  shown 
in  this  section  is  that  any  quantity  which  is  an  average  over  a 
thermal  distribution  can  be  expressed  as  a ratio  of  two  partition 
functions.  To  reiterate  the  general  nature  of  this  result  consid- 
er some  arbitrary  dimensionless  quantity  B,  whose  magnitude  in  the 
state  i is  B^.  The  thermal  average,  <B>  of  B is 

<B>  = ZpiBi.  (2.40) 

It  follows  by  explicit  differentiation  of  (2.40)  that 

RT23<B>/9T  = <BE>-<B><E>  (2.41) 

and  hence  that  an  'activation  energy'  can  always  be  defined 

RT  91n<B>/9T  * -fgj-  - <E>  * <E>g-<E>  (2.42) 

where  <E>0  is  the  average  energy  of  the  normalised  distribution 
bi  = giBiexp(-Ei/RT)/QB  (2.43) 

with 

Qfi  = Zg^expt-E^RT)  . (2.44) 

Recalling  that  pi  = g^exp (-E^/RT) /Q  we  have  from  (2.40)  and  (2.44) 
<B>  = Qb/Q  (2.45) 

The  expression  of  the  reactivity  in  the  form  (2.16)  is  thus 
a mathematical  identity.  Even  if  the  reactants  were  state  select- 
ed, one  could  still  obtain  a suitably  generalised  similar  result. 
The  merit  of  the  observation  is  thus  entirely  in  the  question  of 
whether  it  can  be  used  to  any  advantage  in  extending  our  under- 
standing of  chemical  reactivity. 
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III.  Cycles. 

The  ultimate  test  of  the  concept  of  free  energy  of  the  re- 
active reactants,  Gr=-RTlnQr  , is  whether  it  can  be  employed  to 
construct  thermodynamic  cycles.  The  typical  cycle  to  be  considered 
will  be  a work  cycle  in  which  a system  coupled  to  a heat  bath  is 
restored  to  its  initial  state.  The  second  law  implies  that  such 
a process  cannot  deliver  any  net  work.  We  shall  show  that  the 
definition  of  Gr  is  consistent  with  this  implication  of  the  second 
law,  and  explore  a few  cycles  of  possible  interest.  A general 

demonstration  that  thermodynamic  cycles  can  be  constructed  also 
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for  systems  in  disequilibrium  has  been  given  elsewhere. 

Reactive  Reactants.  The  simplest  possible  cycle  is  one  where  a 
mole  of  reactants  (at  equilibrium  with  a heat  bath)  is  separated 
into  pr  moles  of  reactive  reactants  and  Pnr=l~Pr  moles  of  non- 
reactive reactants.  The  work  of  separation  must  then  be  equal  to 
the  maximal  work  that  can  be  recovered  when  the  two  subsystems  are 
allowed  to  return  to  equilibrium  (in  an  isothermal  process  at  the 
same  temperature  T) . The  proof  is  an  immediate  implication  of  our 
basic  conclusion  (cf.  (2.14))  Q=Qr+Qnr  . Since  pr=Qr/Q  and  since 
the  work  available  from  pr  moles  of  reactive  reactants  (cf.  (2.38)) 
is  Pr(Gr~G)  and  similarly  for  the  non-reactive  reactants,  it  follows 
that 

Wsep=RTI-PrlnPr-pnrlnpnr1-Pr(Gr-G>+Pnr(Gnr-G> • (3'1> 

Nascent  Products.  To  elaborate  on  the  previous  cycle,  let  the 
Pr  moles  of  reactive  reactants  undergo  their  single  collision  pro- 
ducing pr  moles  of  nascent,  unrelaxed,  products.  Work  is  then 
obtained  as  the  nascent  products  are  allowed  to  relax,  isothermally , 


to  equilibrium,  producing  pr  moles  of  thermal  products.  To  con- 
vert these  into  thermal  reactants  requires  -prAG°  units  of  work. 
The  net  work  delivered  by  the  cycle  is 


W = “W  + V*  + Wn  + pAG. 
sep  na  n r 


(3.2) 


Here  = p (G  - G)  is  the  work  available  from  the  p moles 
na  nr  nr  nr 

of  nonreactive  reactants  that  are  brought  back  to  equilibrium  and 

Wn  is  the  work  available  form  the  pr  moles  of  nascent  products. 

To  show  that  W = 0 one  must  show  that  W + p AG°  = W = p (G  - G) 

n r r r r 

where  W is  the  work  available  from  p moles  of  reactive  reactants, 
r r 

for  then,  using  (3.1),  W = 0.  Since  AG°  = (G'  - G)  where  G'  is 

the  free  energy  of  the  thermal  products,  we  need  to  show  that  Gn, 

the  free  energy  of  the  nascent  products  (Wn  = Pr(Gn  - G'))  is  given 


G = G 
n r , 


(3.3) 


i.e.,  the  free  energy  of  the  nascent  products  equals  that  of  the 

reactive  reactants.  Of  course,  one  can  argue  that  (3.3)  must  be 

40 

the  case,  as  otherwise  the  second  law  will  be  violated.  In  other 

words,  if  (3.3)  fails,  one  could  employ  a velocity  selector  (or 

any  other  device  for  selecting  more  reactive  molecules  and  thereby 

affecting  at  least  a partial  separation)  to  produce  work  from  heat, 

on  a continous  basis.  One  can,  however,  prove  (3.3)  without  appeal 

to  Maxwell's  demon,  as  follows: 

To  prove  (3.3)  one  needs  to  invoke  detailed  balance  twice. 

It  is  important  to  note  that  since  we  are  dealing  here  with  isolated 

41 

binary  collisions,  detailed  balance  follows  from  microscopic  re- 
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versibility  and  does  not  depend  on  invoking  the  second  law  in  its 
proof.  The  first  implication  of  detailed  balance  is  that 


Gr  * g; 


(3.4) 


i.e.,  that  the  free  energy  of  the  reactive  reactants  in  the  forward 
reaction  equals  that  of  the  reactive  reactants  in  the  reversed 
reaction.  A simple  short  proof  is  based  on 


(Q'/Q)  = K(T)  = pr/p^. 


(3.5) 


where  K (T)  is  the  equilibrium  constant.  Hence  Qr  = Q^.  and  (3.4) 
follows.  The  second  implication  of  detailed  balance  is  that 


Gn  - g; 


(3.6) 


i.e.,  that  the  free  energy  of  the  nascent  products  of  the  forward 

reaction  equals  the  free  energy  of  the  reactive  reactants  in  the 

reversed  reaction.  The  proof  is  immediate,  since  detailed  balance 
41 

implies  that  the  two  distributions  are  identical.  Combining 
(3.6)  and  (3.4)  we  have  (3.3)  and  hence  conclude  that  W = 0 (cf. 
(3.2))  . 

There  are  altogether  four  distributions,  two  on  either  side 
of  the  barrier.  The  two  on  the  same  side  (reactive  reactants  and 
nascent  products  of  the  same  species  but  of  opposite  reactions) 
are  identical  and  all  four  have  the  same  free  energy.  It  is  this 
last  conclusion  which  enables  one  to  incorporate  kinetic  concepts 
into  thermodynamic-like  cycles. 


I 


Energy  Disposal.  The  free  energy  of  a mole  of  thermal  reactants 
exceeds  that  of  a mole  of  thermal  products  by  G - G'  = -AG°.  On 
the  other  hand,  the  free  energy  of  a mole  of  nascent  products 
exceeds  that  of  a mole  of  thermal  products  by 


G-  G ' =G-G+G-G'  = G-  AG°. 
nr  a 


(3.7) 


If  a reaction  does  not  take  place  spontaneously,  it  follows  that 
GR  = G^,  > G'.  It  takes  work  to  drive  the  reaction  and  the  excess 
of  the  work  required  (G  ) over  the  minimal  work  required  to  form 
the  thermal  products  (AG°)  appears  as  the  work  available  from  the 
nascent  products. 

The  inertia,  G=  - AG°,  is  thus  the  work  available  from  the 

a 

nascent  products  as  these  relax  to  thermal  equilibrium.  Given  two 
reactions  at  the  same  AG°,  the  one  leading  to  more  'work-rich' 
products  will  have  a higher  G and  hence  a lower  reactivity.  For 

a 

a given  AG°,  specificity  of  energy  disposal  as  measured  by  Gn  - G' 

increases  as  the  inertia  increases  and  hence  as  the  reactivity 
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diminishes. 


What  is  still  missing  is  the  dependence  of  the  inertia  on 
AG°.  For  an  exoergic  reaction,  where  -AG°  > 0,  it  follows  from 
(3.7)  that  Gn  * Gr  > G = G'  - AG°,  i.e.,  that  the  excess  free 
energy  of  the  nascent  products  (Gr  - G')  exceeds  -AG°.  However, 
for  an  endothermic  reaction,  where  -AG°  <0,  (3.7)  provides  no 
useful  bound.  What  we  need  is  the  dependence  of  Gn  - G'  on  AG, 
i.e.,  what  is  the  fraction  of  G&  that  ends  up  as  excess  free  energy 


of  the  nascent  products.  Any  cycle  that  we  were  able  to  come  up 
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with,  turned  out  to  be  equivalent  to  (3.7)  and  hence  failed  to 
provide  an  answer.  The  question  will  thus  have  to  wait  till  sec- 
tion VI  where  the  conclusion  will  be  that  G - G'  - -AG°  (i.e., 

n 

Ga  * 0)  as  -AG°  >>  RT  while  G^  - G'  =0  (i.e.,  G&  = AG^  as  AG° » RT. 
Reactivity  as  an  Equilibrium  Constant.  The  introduction  of  cycles 
enables  one  to  regard  Pr(T)  = Q^/Q  as  the  equilibrium  constant 
between  the  thermal  reagents  and  the  reactive  reagents,  say 

Cl  + HI  £ (Cl  + HI) r . 

As  a check  we  verify  that  it  satisfies  the  usual  product  rule  for 
a sequence  of  equilibria.  Consider  say 

Cl  + HI?=iClH  + I. 

It  can  be  broken  into  the  following  sequential  stages: 

Cl  + HI^=*(C1  + HI)r^=5(ClH  + I)n5==*ClH  + I. 

The  respective  equilibrium  constants  are  pr(T),  unity  (since  Gr  = 
Gn)  and  l/p^(T) . Indeed  since  Qn  = Qr  , 

K(T)  = pr(T).l.(l/p^.(T))  . (3.8) 

As  another  consistency  check,  it  follows  that  reactive  reactants 
have  a unit  reactivity,  as  they  should. 

Branching  Ratios.  Consider  first  an  intramolecular  branching  ratio 

43 

as,  say,  for  the  F + HD  reaction 
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F + HD 


FH  + D 


\ 


FD  + H 


The  relative  rate  of  H vs.  D abstraction,  starting  with  thermal 
reagents,  is  given  by 


-RTlnr  = -RTln(p  /p  ) - G (H)  - G (D) 
h u n n 


(3.9) 


where  Gn(X)  is  the  free  energy  of  the  nascent  products  in  the  X 


abstraction  path.  T is  thus  the  ratio  of  the  partition  function 


of  the  nascent  products  (or  of  the  reactive  reactants) , T = Q (H)/Q  (D) 

n n 


Alternatively,  if  Kn(T)  is  the  equilibrium  constant  between  F + HD 


and  say  (FD  + H)n  (i.e.,  Kn(T)  = PD(T))  then 


log  T = APKn* 


(3.10) 


An  explicit  cycle  corresponding  to  the  identity  (3.10)  is  easily 
constructed.  Consider  the  two  reactions 


, D + HF  H + DF 

H + DFT  and  D + HF^ 

kHD  + F ^F  + HD 


whose  branching  ratios  (for  thermal  reagents)  will  be  denoted  by 


rR  and  r respectively.  Then  (3.10)  is  equivalent  to 


r * VrD  * 


(3.11) 
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The  identity  (3.10)  remains  equally  valid  also  for  comparisons. 

r 

A very  simple  example  is  a comparison  of  the  rates  of  forward  and 
reversed  reactions.  Using  the  cycle 


A + BC  -►  (AB  + C) 


(A  + BC)n  ■*-  AB  + C 


and  (3.4),  we  recover  the  well  known  result  kr/k^.  = K.  More  elab- 
orate cycles  can  also  be  devised.  An  alternative  direction  is, 
as  in  the  intramolecular  branching  problem,43  to  relate  the  dif- 
ferences in  reactivity  to  the  energy  disposal  (or  energy  consump- 
tion) . Two  worthwhile  examples  are44 

Cl  + HI  -*■  C1H  + I vs.  Cl  + HBr  -►  ClH  + Br 


H + 02N  + HO  + ON  vs.  H + 03  HO  + 02  . 


For  both  examples,  the  first  reaction  is  over  an  order  of  magnitude 
more  reactive  than  the  second.  In  section  VI,  where  the  dependence 
of  G&  on  AG°  is  explored,  it  will  be  found  possible  to  obtain  ex- 
plicit results  (cf.  figure  3 below)  for  the  intramolecular  branching 
ratio. 


IV.  The  Partition  Function 


A direct  computation  of  the  reactivity  requires  the  determin- 
ation of  the  partition  function  Qr.  To  do  so  we  need  the  'degenera- 
cies' gri  of  the  reactive  reactants  (cf.  (2.22))  and  these  are  only 
available  from  a detailed  microscopic  (experimental  or  computational) 
study.  This  section  explores  two  theoretical  techniques  that  may 

prove  useful  for  determining  Qr  in  terms  of  macroscopic  observables . 
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The  first  is  the  use  of  the  Legendre  transform.  This  venerable 
procedure  will  be  invoked  below  in  the  discussion  of  free  energy 
relations.  The  second  is  the  use  of  bounds. 

The  Legendre  Transform.  The  partition  function  is  a function  of 

the  temperature.  The  Legendre  transform  replaces  it  by  a different 

function,  whose  variable  is  -31nQr/3 (1/T)  such  that  both  functions 
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provide  precisely  the  same  information.  Using  (2.33)  the  Leg- 
endre transform  of  lnQr  with  respect  1/T  is  Sr 

sr  = lnQr  - (1/T) [-31nQr/3(l/T)]  (4.1) 

Knowing  Sr  as  a function  of  <E>r  is  thus  equivalent  to  knowing  Qr 
as  a function  of  T.  Indeed,  the  determination  of  Sr  by  integrating 
(2.36)  was  essentially  the  suggestion  of  LaMer.®  Unfortunately, 


the  determination  of  Qr  via  the  measured  T-dependence  of  the 
activation  energy  is  not  easy  to  implement,  due  to  experimental 
difficulties  in  detecting  small  changes  in  E . Similar  practical 
difficulties  arise  in  another  attempt  to  invert  a Legendre  trans- 
form of  Qr,  to  be  discussed  below.  Much  earlier,  such  difficulties 
in  the  determination  of  Q have  given  rise  to  the  formulation  of 
the  Nernst  heat  theorem.  There  may  well  be  a lesson  here  too,  but 
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its  significance  is  not  clear  at  the  moment.  ' 


Energy  Consumption.  Qr  may  well  depend  on  additional  variables 
beside  the  temperature  via  its  dependence  on  the  gr^’s.  The  ratio 
gri^9i  *s'  construction,  the  probability  of  reaction  for  a state 
in  the  group  i and  the  observation  that  gr^  is  not  necessarily 
proportional  to  g.  (i.e.  g^/g^  7*  Pr)  gives  rise  to  the  concept  of 
’selectivity  of  energy  consumption’.  To  allow  also  for  the  possib- 
ility that  the  probability  of  reaction  depends  not  only  on  the 
total  energy  but  also  on  the  partition  of  the  total  into  the  differ- 
ent degrees  of  freedom,  we  shall  regard  i as  an  index  of  a detailed 
state  selection  for  the  reactions.  To  bring  this  aspect  into  the 
open  it  is  convenient  to  introduce  S [x] , 

(4.2) 

the  entropy  of  the  distribution  of  the  reactive  reactants,  regarded 


S [x]  = -Ex  . ln(x  ./g. ) , 


S [x]  = S -Ex  .ln(g  ,/q.) 

r i n ^ri  ’i 


as  a distribution  of  quantum  states.  S [x]  is  maximal  when  all 

4 

states  at  a given  energy  react  with  the  same  probability.  From 
their  respective  definitions, 

(4.3) 

and  since  gri  < g^,-ln (gri/g^) £0  with  equality  if  pf  = 1 and  so 

S[x]  Sr  with  equality  when  Pr  = 1.  It  follows  that 

pr  < exp{S [x]-<E>r/RT}/Q.  (4.4) 

Two  reactions,  with  the  same  energy  of  activation,  may  thus  differ 

in  their  reactivity  because  of  differences  in  their  selectivity  of 

energy  consumption.  As  the  selectivity  increases,  S [x]  diminishes 
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and  so  does  the  reactivity.  The  origin  of  this  ’selectivity- 
reactivity’  principle  is  indeed  trivial.  Since  at  most  all  the 


trajectories  in  the  group  can  be  reactive,  gr^  < gr . Since 
selectivity  implies  that  gr^  ^ prg^  follows  that  a measure  of 
selectivity  should  vanish  when  g^  = and  otherwise  be  always 

positive.  Using  (2.12),  the  'local'  deviance  of  gri  from  Pr9^  is 
(gri/Pr9i) = (Xri/Pi) • ^ 'integral'  measure  of  selectivity  is  thus 

given  by  qr  defined  as 

lnq  = -Ex  .ln(x  ./P-)  (4.5) 

^r  . n n ri 

1 4 

where,  using  the  Gibbs-Boltzmann  inequality  , lnqr  < 0 with  equality 

if  and  only  if  gri  = gjPr.  However,  using  (2.22) 

lnq  = S [x] - ( <E>  /RT) -InQ  (4.6) 

r r 

or 

q^  = exp{S [x]-<E>r/RT}/Q  (4.7) 

so  that  pr  < qr . This  bound  on  the  reactivity  will  be  further 
discussed  below. 

A Tale  of  Two  Entropies.  For  a molecular  system,  entropy  is  a 

measure  of  the  uncertainty  as  to  the  actual  quantum  state  of  the 
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system  . The  magnitude  of  the  entropy  will  thus  depend  on  the 
amount  of  information  already  available  on  the  system.  Any  input 
that  can  be  used  to  eliminate  some  quantum  states  as  possible 
states  will  thus  tend  to  reduce  the  entropy.  The  inequality 
S [x]  » Sr  is  a quantitative  statement  of  such  consideration  S [x] 
is  the  uncertainty  given  that  the  distribution  of  quantum  states 
is  xri-  The  choice  is  then  among  all  possible  states  of  the 
reactants.  The  work  required  for  such  a search  is29,30  RTlnqr. 

Sr  is  the  uncertainty  given  that  the  distribution  is  xri  and  that 
the  choice  is  restricted  to  the  reactive  reactants.  Since 
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gri  < g^f  the  number  of  possible  quantum  states  is  reduced, 

Sj,  « S[x]  and  the  work  required  for  such  a search,  RTlnpr»  is  lower. 
The  difference  between  the  two  entropies  is,  in  a sense,  a quantal 
effect.  In  the  classical  limit  we  can  take  i to  index  the  initial 
conditions,  whence  gr^  is  either  zero  or  unity,  while  g^  is  unity. 
When  gri  = 0 xri  also  vanishes, so  that  the  two  entropies  agree,  as 
the  search  is  precisely  over  the  same  range  of  possible  contenders. 
The  same  conclusion  is  quantitatively  evident  from  (4.3).  Of  course, 
the  need  for  two  entropies  is  partly  of  our  own  making,  for  having 
introduced  the  distribution  of  reactive  reactants.  In  general, 
it  is  S[x]  which  has  the  conventional  structure  of  an  entropy,  for 
it  is  the  missing  information  of  the  distribution  xr^,  given  only 
macroscopic  information  on  that  distribution.  It  is  therefore  of 
interest  to  note  that  the  can  be  obtained  as  a Legendre  trans- 
form S [x]  . 

Consider  the  simplest  situation  where  besides  <E>r»  the  only 
macroscopic  input  required  to  characterise  the  distribution  of 
reactive  reactants  is  the  magnitude  of  <C>r  = Ix^C^  where  is  a 
property  which  depends  on  i.  Then 

xri  = gi  exp(- (Ei/RT)-YCil/Qr  (4.8) 

where  the  magnitude  of  y is  adjusted  so  that  the  distribution  xr^ 
is  consistent  with  the  value  of  <C>r.  For  that  value  of  y, 

<C>r  = -31nQr/3y.  (4.9) 

Since  gri  = g^exp(-yC^),  it  follows  that  yC^  > 0 and  y<C>r>0  with 
equality  if  gri  = gi. 

The  Legendre  transform  of  lnQr  with  respect  1/T  is  still 
Sr,  but  the  Legendre  transform  with  respect  both  1/T  and  y is 
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S [x]  = lnQ-(l/T)31nQr/3(l/T)-Y31nQr/3Y  (4.10) 

= InQ  + <E>r/RT  + Y<C>r 
= Sr  + Y<C>r 

Here  S[x]  is  a function  of  <E>r  and  <C>r#  whereas  Sr  is  a function 
of  <E>r  and  of  y.  Since  y - 3S[x]/3<C>r,  Sr  and  S [x]  are  also  a pair 


of  Legendre  transforms 

Sr  - S[x]-<c>  ||isi 


(4.11) 


Bounds . The  upper  bound  (4.4)  on  the  reactivity  was  derived  from 
the  observation  that  not  all  trajectories  are  necessarily  reactive. 
This  suggests  a generalisation.  There  is  a range  of  intermediate 
distributions  between  p^  ('all  trajectories  that  originate  from  the 
reactants')  and  x ^ ('reactive  trajectories').  For  example,  all 
trajectories  that  originate  from  the  reactants  and  that  cross  any 
dividing  surface  which  separates  the  reagents  from  the  products.47-50 
If  the  surface  is  close  to  the  reactants,  the  distribution  is  close 
to  p^  while  if  it  is  close  to  the  products,  it  is  close  to  xr^. 

If  there  is  a 'surface  of  no  return'  such  that  any  trajectory  which 
crosses  the  surface  proceeds  to  form  products,  then  the  distribution 
for  that  surface  is  precisely  x In  mathematical  terms,  we  argue 


that  since  gri  « g^  one  can  always  define  g+^  such  that 

g . < g, . < g.  (4.12) 

’ri  ^+i 

Multiplying  every  term  in  (4.12)  by  the  (positive)  number  exp (-E^/ FT) 
and  summing  over  i, 

Qr  < Q+  « Q (4.13) 


Pr  < p+  = Q+/Q 


(4.14) 
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The  physical  argument  is,  of  course,  the  one  employed  in  the 

47-49 

variational  theory  of  reaction  rates  , and  enables  one  to  endow 
the  mathematical  result  (4.14)  with  a physical  interpretation. 

A mathematical  proof  of  the  bound  is  also  available  using  the 
Gibbs-Boltzmann  inequality  as  in  (4.4). 

A lower  bound  for  the  partition  function  is  also  readily 
derived,  by  a simple  interchange  of  xr^  and  p^  in  the  Gibbs-Boltzmann 
inequality 

Zpiln(pi/xri)  = -S  [p] - ZpJ.lnCxri/gi)  (4.15) 

= -S[p]+<E>/RT  + y<C>  + lnQr 
= lnQr  - InQ  + y<C>  > 0 

or  lnpr  > -y<C>.  The  two  inequalities  brought  together  are 

exp[-Y<C>]  « Pr  « exp  [~Y<C>r]  (4.16) 

where  <C>  is  the  average  value  of  (cf.  (4.8))  over  the  distri- 
bution of  thermal  reactants. 

The  link  between  the  upper  bound  in  (4.14)  and  (4.16)  is 
provided  by  noting  that  the  essential  inequality  is  Qr<Qexp [-Y<C>r] 
and  that  using  (4.12)  one  can  also  write  Qr  < Q+exp [~y (<C>r“<C+>r) ] 
where  -ln(g+i/gi)  « YC^,  so  that  Q . « Q+  and  especially  obtains 

when  there  is  a surface  of  no  return  so  that  ct=C^.  A similar  trick 
is  available  for  the  lower  bound,  Qr  > Q+exp [~y (<C>+-<C+>+) ] where 
the  + subscript  denotes  an  average  over  the  distribution  g+^  exp 
(-Ei/RT)/Q+.  Since  Y<C>+5>  0,  one  can  obtain  a lower  bound 

Pr  » (Q+/Q)  exp (y<C+>+)  (4.17) 

which  is  complementary  to  the  upper  bound  (4.14).  Since  yC^  = 
-ln(g+^/gi),  the  bound  (4.17)  can  be  computed  when  Q+  is  available. 
The  bound  is  not  perfect,  since  it  is  always  a bound  rather  than 


t 
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an  equality,  but  it  can  be  improved  upon. 
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V.  Dynamics 

This  section  introduces  the  concept  of  reactivity  on  collision- 
theoretic  grounds.  As  in  section  II,  a classical  trajectory  compu- 
tation point  of  view  is  employed. 

Let  n denote  the  set  of  classical  variables  corresponding  to 
the  quantal  interval  quantum  numbers,  the  translational  energy 
of  the  collision  and  b the  impact  parameter,  whose  range  is  taken 
from  zero  to  some  maximal  value  B.  The  number  of  trajectories 
used  in  the  range  b to  b + db  is  related  to  the  total  number  run  by 

N(n,ET,b)db  = (2irbdb/irB2) N (n,ET>  . (5.1) 

Here  N(n,ET)dET  the  number  of  trajectories  used  for  translational 

energy  in  the  range  ET,  E^  + dE^  with  specified  n.  The  number  of 

reactive  trajectories  for  an  impact  parameter  in  the  range  b,  b + 

33 

db  is  Nr (n,ET,b)db.  The  reactivity  at  a given  b (sometimes  called 
'the  opacity  function')  is 


P(b)  = Nr(n,ETb)/N(n,ET,b) , (5.2) 

and  the  number  of  reactive  trajectories  for  ET  in  the  range  E^, 

Et  + dET  and  given  n is  a sum  (i.e.,  integral)  over  all  possible 
values  of  b 


-27- 


The  reaction  cross  section,  defined  by33 
B 

a (n,E  ) =2 it  f bP(b)db  (5.4) 

T o 

can  now  be  rewritten,  using  (5.2)  and  (5.1)  as5^ 

(5.5) 

? B 2 
a(n,ET)  = f Nr(n,ET,b)db/N(n,ET)  = i:B^Nr  (n,ET) /N  (n,ET) 

o 

2 

i.e.,  as  irB  times  the  reactivity  at  a given  ET.  A similar  'Monte- 
Carlo'  integration  can  also  be  used  to  relate  the  rate  constant 
and  the  reactivity  at  the  temperature  T.  To  do  so  we  sample  tra- 
jectories at  different  translational  energies  by 

N(n,ET)  = Z(ET)f(ET)N(n)/Z  (5.6) 

Here  N(n)  is  the  number  of  trajectories  employed  for  reactants 
in  the  internal  state  n.  f (E^) dE^  is  the  fraction  of  systems  at 
thermal  equilibrium  with  translational  energy  in  the  range 
Et  ,Et  + dET 

(5.7) 

f(ET)  = Pt(Et)  exp(-ET/RT)/Q^ 

where  P^(E^,)  is  the  translational  density  of  states  and  Qfc  is  the 

translational  partition  function.  Z(E^)  is  the  collision  frequency 

2 

at  the  translational  energy  E^,,  Z(ET)  = uB  v,  where  v is  the  relative 

velocity  and  Z is  the  thermal  average  of  Z(E,p).  If  B is  independent 
2 

of  E_,  Z = ttB  <v>. 

T — 

The  reaction  rate  constant  is 
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k(n,T)  = f f (ET) va(n,ET)dET 


(5.8) 


and  using  (5.5)  and  (5.6), 

00 

k(n,T)  = / dETf(ET)Z(ET) [Nr(n,ET)/N(n,ET)]  = Z [N  (n) /N (n) ] . 

o 

(5.9) 

Here  Nr(n)  is  the  number  of  reactive  trajectories 


Nr(n)  = / dETNr(n,ET), 


(5.10) 


starting  with  the  internal  state  n. 

The  final  stage  is  to  obtain  the  thermal  rate  constant 


k(T)  = E [N(n)/N]k(n,T)  = Z (N  /N) 


(5.11) 


n 


where  N = E N(n)  is  the  total  number  of  reactive  trajectories, 
r n r 

The  careful  reader  will  note  that  Z,  the  collision  rate,  de- 
pends on  B,  whose  numerical  magnitude  has  not  been  precisely  assigned. 
This  indeed  is  a weakness  of  classical  mechanics.  It  is  unable  to 
specify  a B beyond  which  there  are  no  collisions.  Even  so,  the 
magnitude  of  k(T)  is  independent  of  the  choice  made  for  B (provided 

only  it  is  large  enough,  and  that  N is  large) . The  reason  is  that 

2 

as  can  be  inferred  from  (5.1),  it  is  only  the  ratio  B /N  that  matters 
and  only  this  ratio  appears  in  k(T) . 


The  Sampling  Procedure.  It  is  essential  to  note  that  by  running 
enough  trajectories  (i.e.,  in  the  limit  of  large  N)  one  obtains 
the  correct  (quasi-classical)  magnitude  of  the  rate  constant  for 
any  reasonable  sampling  strategy  of  the  initial  states.  The  reason 
is  obvious.  If  we  accentuate  the  role  of  any  particular  state, 
the  number  of  reactive  trajectories  originating  from  that  state 
will  also  increase,  and  the  reactivity  (cf.  (5.2)  or  (5.5))  will 
remain  unchanged,  provided  N is  large.  Thus,  while  the  functional 
form  for  the  rate  constant  (as  the  collision  frequency  times  the 
fraction  of  trajectories  that  do  react)  does  depend  on  a particular 
choice  of  the  sampling  procedure,  its  magnitude  computed  using 
(5.8)  is  independent  of  this  choice.  In  other  words,  the  ratio 
Nr(n,E)/N(n,E)  (which  is  the  analogue  of  - gri^gi  of  sec“ 

tion  II)  is  independent  of  the  sampling  procedure.  A different 
sampling  procedure  is  employed  in  the  Appendix  to  rederive  the 
present  results  and  to  bring  them  to  the  form 

k(T)  = <(T)kTST(T)  (5.12) 

Here  kTST(T)  is  the  transition  state  theory  rate  constant  and  K(T) 

< (T)  = [Nr/N+]  < 1 (5.13) 

is  the  transmission  coefficient.  N+  is  the  number  of  trajectories 
that  originate  from  the  reagents  and  that  cross  the  dividing  surface 
from  reagents  to  products  and  hence  is  an  upper  bound  for  Nr»  As 
is  already  evident  and  will  be  further  discussed  in  the  appendix, 
k(T)  can  be  regarded  as  the  reactivity  of  the  transition  state. 
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VI.  The  Brtfnsted  Series 

The  variation  of  the  magnitude  of  the  reaction  rate  constant 
with  structural  modifications  in  the  reagents  is,  in  principle, 
contained  in  the  two  partition  functions  Qr  and  Q.  In  practice, 
equilibrium  statistical  mechanics  accounts  for  the  variations  in 
Q,  but  we  have  no  obvious  route  (except  for  cycles)  to  Qf.  A simple 
approach,  based  on  the  concept  of  a Br^nsted  series,  is  examined 
in  this  section,  where  the  point  of  view  is  strictly  operational. 

An  interpretation  of  the  results  in  terms  of  a model  is  considered 
in  the  next  section. 

The  BrgSnsted  Coefficient.  The  Br^nsted  coefficient  is  a measure 
of  the  change  in  the  reaction  rate  constant  due  to  changes  in  the 
free  energy  of  the  reaction 


a = 31nkr(T)/31nK(T)  = -31nkr (T) /3 (AG°/RT) . (6.1) 


The  mathematical  definition  is  given  a physical  interpretation  by 

the  assumption  that  in  a series  of  'similar'  reactions,  structural 

changes  in  the  reactants  change  the  reaction  rate  primarily  via 

the  change  in  the  AG°  of  the  process.  This  interpretation  is  saved 

from  being  a circular  one  ('a  series  of  similar  reactions  is  one 

for  which  the  Br^nsted  relation  is  a useful  measure',  'the  Br^hsted 

slope  is  a useful  measure  for  a series  of  similar  reactions')  by 

the  empirical  observation  that  there  are  indeed  many  known  examples 

2 18  28 

of  the  utility  of  the  concept.  ' ' The  very  same  observation 

is  valid  of  course  in  general  ('The  phase  rule  is  valid  for  a hetero- 
genous system  in  equilibrium',  'a  heterogenous  system  is  in  equili- 


brium if  the  phase  rule  obtains').  The  only  logical  grounds  for 


-31- 

further  enquiry  are  thus  that  one  would  like  to  link  the  definition 
to  other,  possibly  more  general,  physico-chemical  concepts  and, 
in  so  doing,  identify  the  implications  of  deviant  behavior.  In 
other  words,  one  would  like  to  invoke  the  concept  as  a yardstick 
in  the  same  fashion  that  the  phase  rule  is  used  to  decide  whether 
geological  samples  (or  meteorites)  are  indeed  in  equilibrium. 

The  'Entropy'  Term.  For  small  variations  in  AG°,  the  variation 
in  the  rate  constant  can  be  expressed  as11 

61nkr  = -ct6(AG°/RT)  (6.2) 

However,  if  a is  not  independent  of  AG°,  the  application  of  (6.2) 
is  limited  (say  to  changes  in  AG°  due  to  external  parameters,  e.g. , 
pressure).  In  particular,  structural  modifications  such  that  the 
change  in  AG°  exceeds  RT  cannot  be  discussed.  An  alternative  ap- 
proach is  to  explicitly  integrate  (6.1),  holding  all  variables 
except  AG°  constant.  To  allow  for  a possible  AG°  dependence  of 
a,  (6.2)  is  written  as 

dlnkr(T)  = -d (aAG°/RT)  + (AG°/RT) da  (6.3) 

or 

Ink  (T)  = InZ  - a (AG°/RT)  - (G°/RT)M(a)  (6.4) 

jl  a 

Here  InZ  is  the  integration  constant , which  is  thus  independent 
of  a while  M(a)  is  a dimensionless  function  of  a,  defined  by  com- 
paring (6.4)  to  (6.3),  as 

(6.5) 

rA 


(G°/RT) 9M(a)/3a  = -AG°/RT. 

a 


f 
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In  other  words,  G°M(a)  is  defined  as  the  integral  of  AG°dot.  G 
has  the  dimension  of  free  energy  and,  by  construction,  is  indepen- 
dent of  a.  A Brtfnsted  series  is  characterized  by  a given  value 
of  G°.  For  such  a series  one  can  explicitly  construct  the  function 

cl 

M(a) , given  the  magnitude  of  the  reaction  rate  for  different  members 
of  the  series.  (This  serves  to  determine  a as  a function  of  AG°, 
and  equation  (6.5)  can  then  be  explicitly  integrated.)  The  magnitude 
of  G°  can  be  fixed,  say  by  the  convention  that  at  the  stationary 

8i 

point  of  M(ct),  (that  is,  at  G°3M/3a=-AG°  = 0),  G°  = G . 

a.  a a 

The  representation  of  the  reaction  rate  constant  as 

k (T)  = Z exp (-G  /RT)  (6.6) 

a 

(6.7) 

with  Z independent  of  AG°  is  thus  a mathematical  statement,  equiva- 
lent to  the  original  definition  of  a.  Given  equation  (6.2),  there 

always  exists  a function  G°M(a)  such  that  the  AG°-dependence  of 

a 

the  rate  constant  can  be  stated  in  the  form  (6.6)  with  G=  having 
the  form  (6.7),  and  using  the  definition  (6.5),  satisfying  the 
Legendre  transform  relations 

(G  /G°)  = M(a)  - a (3M/3a) 
a a 

or 

M(a)  = (G  /G°)  - AG°[3  (G  /G°)/3AG°) 

a a ci  cL 


The  specification  of  the  a-dependence  of  M is  thus  strictly 

equivalent  to  the  determination  of  G as  a function  of  AG°.  (An 

a 


(6.8) 

(6.9) 


G = oAG°  + G°M(a) 
a a 
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explicit  example  is  discussed  later.)  For  small  variations  in 
AG°,  the  result  (6.2)  remains  valid  for  any  function  M(a)  that 
satisfies  (6.5)  or,  equivalently,  (6.8),  but  only  if  G°  >>  AG° 

cL 

can  a be  taken  to  be  independent  of  G. 

Detailed  Balance.  To  inject  some  physics  into  the  discussion, 
consider  the  implication  of  detailed  balance.  For  a given  ele- 
mentary reaction,  the  forward  and  reverse  rate  constants  are  re- 
lated by 

k(T)/k'(T)  = K (T)  = exp(-AG°/RT) . (6.10) 

Since  Z (cf.  equation  (6.6))  is  independent  of  AG°,  it  can  be  taken 
to  be  the  same  for  the  forward  and  reverse  reaction  (always  subject 
to  the  assumption  that  these  belong  to  the  same  Br^nsted  series) . 

Hence  the  Br^nsted  slope  for  the  reversed  reaction,  (whose  standard 
free  energy  change  is  -AG°,  where  AG°  is  the  standard  free  energy 
change  in  the  forward  reaction) , is 

a'  = -91nk ' (T) /9 (-AG°/RT)  = -91nK (T) /9 (AG°/RT)  + 9 Ink (T) /9 (AG°/RT) 
= 1 - a 
or 

a + a'  = 1.  (6.12) 

Furthermore,  since  (6.10)  implies  that 


G (a)  - G (a' ) = AG° 

a a 


(6.13) 


1 


I 

M(a)  - M(a')  = 0,  or  using  (6.12), 

M(a)  = M ( 1 - a) . (6.14) 

As  a varies  over  the  range  zero  to  one,  M(a)  is  a symmetric  func- 
tion of  its  argument.  Indeed,  we  know  from  (6.5)  that  M(a)  has 
a stationary  point  (9M(a)/3a  = 0)  for  that  value  of  a that  cor- 
responds to  AG°  = 0 (a  'symmetric'  reaction).  The  symmetry  of 
M(a)  implies  that  the  maximum  occurs  at  a = 1/2.  On  the  other 
hand,  taking  the  derivative  of  both  sides  of  (6.5)  with  respect 
to  AG°/RT,  we  conclude  that 

(G°/RT) [32M(a)/9a2] [9a/9 (AG°/RT) ] =-l.  (6.15)  j 


The  two  derivatives  in  (6.15)  are  thus  necessarily  of  opposite 
sign.  The  result  is  again  very  much  expected  since  it  is  equi- 
valent to  the  statement  that  the  second  derivatives  of  M(a)  and 
of  G (AG°) , which  are  a Legendre  transform  pair,  are  of  opposite 

a 

signs.  Hence,  if  M(a)  has  a maximum  at  a = 1/2,  a is  a raontonically 
increasing  function  of  AG°. 

The  monotonic  variation  of  a with  AG°  holds  for  any  Br0nsted 
series.  Any  failure  of  this  result  is  thus  evidence  that  one  is 
not  comparing  reactions  of  the  same  series.  All  that  remains  is  to 
determine  whether  a increases  or  decreases  as  AG°  is  increased. 

The  Inertia  of  Chemical  Reactions.  This  is  about  as  far  as  one 
can  go  on  general  grounds.  A small  infusion  of  chemical  common 
sense  (or  empiricism)  is  however  quite  beneficial  at  this  point. 

The  essential  idea  is  due  to  Evans  and  Polanyi^ ' ^ , who  pointed 

A 
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out  that  the  inertia  of  chemical  reactions  exceeds  their  thermo- 
dynamic driving  force.  In  mathematical  terms,  this  means  that  at 
any  magnitude  of  AG°  the  function  G ( AG°)  is  everywhere  above  its 

d 

tangent  at  that  point,  i.e.,  G (AG°)  is  a concave  function  of  AG 

a 

as  is  evident  in  figure  1.  In  more  direct  terms,  consider  a 
symmetric  (AG°  = 0)  reaction  and  another  reaction  of  the  same  family 
with  a small  but  finite  AG°.  Then,  if  (AG°/RT)  > 0,  G increases 

cl 

in  excess  of  AG°/RT  and  if  (AG°/RT)  < 0,  G decreases  by  less  than 

d 

-AG°/RT . These  two  equivalent  statements  are  that  32G  /3(AG°)2  > 0 

cL 

2 2 

and  that  3 M/3a  < 0,  respectively.  Either  one  implies  that 
3a/3AG°  > 0.  Explicitly,  from  (6.15), 

3a/3AG°  = -l/G°[32M/3a2]  > 0.  (6.16) 

d 

Since  a = 1/2  for  AG°  = 0,  it  follows  that  when  AG°  > 0,  a < 1/2 
and  G is  only  moderately  AG°  dependent  whereas  when  AG°  > 0,  G 

cl  cl 

is  both  larger  and  more  strongly  dependent  (a  > 1/2)  on  AG°. 

The  Reactivity-Selectivity  Principle.  The  correlation  of  diminished 
reactivity  (higher  G ) with  increased  selectivity  with  respect  to 

d 

structural  changes  (i.e.,  higher  a = 3G  /3AG°) , is  sometimes  known 

d 
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Example . Figure  3 shows  a plot  of  G&  vs.  AG°  for  the  series  of 

2 8 

gas  phase  reactions  X + H2  -*•  XH  + H (squares)  . Included  in  the 

52 

figure  are  results  (triangles)  for  the  reversed  reactions 
XH  + H -*•  X + H2  obtained  using  (6.10).  X is  a halogen  atom.  The 
continous  curve  is  a fit  to  the  data,  and  a vs.  AG°,  determined 
from  the  continuous  curve,  is  shown  as  the  right  hand  side  pare! 
of  the  figure.  The  Legendre  transform,  M(a)  of  the  continuous 
curve  is 


M(a)  = -[aloga  + (1  - a)log(l  - a) ]/log2,  (6.17) 

so  that 

a = 1/[1  + exp(-AG/G°)  ] , (6.18) 

cl 


is  the  functional  form  shown.  Many  additional  examples  of  the 
use  of  the  functional  form  (6.17)  to  fit  experimental  data  will 
be  found  in  references  25  and  28.  The  results  shown  in  figure  3 
are  however  exceptional  in  the  range  of  AG°  which  is  included. 

For  most  examples,  the  span  in  AG°/G°  is  much  more  limited,  thereby 

d 

preventing  a critical  examination  cf  the  a-dependence  of  M(a) . 

The  problem  is  that  very  general  considerations  suffice  to  show 
that  M(a)  is  symmetric  about  a = 1/2,  where  it  is  maximal.  While 
many  functions  qualify  (e . g. , ' 27 ' 53 • 54  M(a)  = 1 - 4(a  - 1/2)  2 
= 4a (1  - a)),  they  all  must  have  the  same  behavior  about  a = 1/2 
and  hence  can  only  really  differ  in  their  wings,  i.e.,  for  AG° 
substantially  different  from  RT.  Of  course,  this  is  precisely  the 
region  which  is  not  easily  measured  in  the  gas  phase  and  where 
additional  complications  enter  in  solution  kinetics. 
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The  selectivity-reactivity  principle  is  very  much  evident 

in  the  left  panel  of  figure  3.  It  is  worthwhile  however  to  draw 

attention  to  the  scale  on  the  ordinate.  It  is  a logarithmic  scale, 

which  tends  to  gloss  over  small  variations.  Isotope  effects  (i.e., 

replacing  H by  D)  where  the  changes  in  AG°  are  small  cannot  be 

25  2 8 

accounted  for  in  this  fashion,  ' and  do  require  a change  in 
the  magnitude  of  G°.  Similarly,  the  temperature  dependence  of 

a 

the  rate  constant  requires  that  G°  does  depend  on  T.  In  the  more 
general  case  of  reactant  state  selection,  G°  must  be  a function 

cl 

of  the  energy  of  the  reactants,  invoking  detailed  balance,  this 
implies  that  rate  constants  into  different  products  energy  states 
will  have  different  G°  values  and  should  not  be  put  on  a common 

cl 

Br^nsted  plot.  The  discussion  in  section  III  specifies  what  needs 
to  be  done.  M(a)  is  the  Legendre  transform  of  G only  with  respect 
to  aG°.  If  the  role  of  additional  variables  is  to  be  explored, 
one  needs  to  carry  out  a Legendre  transformation  with  respect  to 
these  other  variables  as  well.  While  the  principle  is  thus  clear, 
considerable  work  is  required  in  practice  before  the  utility  of 
such  an  approach  is  established. 

The  Activation  Energy.  A correlation  between  the  activation  energy 
and  the  standard  energy  change  AE°  of  the  reaction  obtains  for  any 
Br^nsted  series.  Applying  the  Gibbs-Helmholtz  equation  (2.31)  to 
(6.7)  and  using  (6.5)  one  finds 


E = -T28 (G  /T)/3T  = aAE°  + E°M(a)  (6.19) 

a a a 

where 

E°  = -T23(G°/T)/3T. 

a SL 


(6.20) 
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In  particular,  it  follows  that  for  a Brjrfnsted  series,  the  height 
of  the  potential  energy  barrier  to  reaction  would  be  correlated 
with  AE°  , the  change  in  bond  energy.  Figure  4 shows  such  a cor- 
relation for  the  same  series  shown  in  figure  3. 

Summary.  A critical  evaluation  of  the  scope  and  limitation  of 
the  functional  representation  G = aAG°  + G°M(a)  requires  that 
the  rate  constants  be  known  for  a series  of  reactions  spanning 
a wide  range  in  AG°/RT.  In  addition,  it  would  be  very  useful  to 
know  these  rates  over  a wide  temperature  range  and  to  have  some 
indication  of  the  dependence  of  the  rates  on  selective  reagent 
excitation  (or  for  different  products  energy  states) . The  two 
series  which  appear  as  likely  candidates  for  such  an  evaluation 

are  X + H2  -*•  XH  + H and  X + HY  -*■  XH  + Y,  where  X and  Y are  halogen 

44  55 

atoms.  What  is  required  is  more  data  on  the  temperature  de- 

pendence of  the  rate  constants.  Two  other  families  which  will 
help  shed  light  on  G°  are56  H + XY  -*■  HX  + Y and57  H + XH  -+■  HX  + H. 

C l 


VII  The  Interpolation  Model 

A simplistic  model  for  the  form  of  the  free  energy  of  activation 

is  introduced  in  this  section.  It  offers  a complementary  point  of 

view  to  the  operational  description  of  section  V I . There  we  took 

it  for  granted  that  a Br^nsted  series  of  reactions  does  exist  and 

concluded  that  this  implies  the  existence  of  a function  M(a)  whose 

properties  were  examined.  Here  we  postulate  an  interpolation 

formula  between  the  free  energies  of  the  reagents  and  the  products 

and  conclude  that  a Br^nsted  series  can  be  defined  and  that  a can 

be  given  a structural  interpretation,  as  the  location  of  the 

12 

’transition  state’  enroute  from  reagents  to  products. 

The  Progress  V ariable.  Consider  n to  be  a progress  variable 

describing  a concerted  chemical  reaction.  n=0  at  the  reagents 

side  and  =1  at  the  products.  It  is  suggestive  to  think  of  n as 
58 

the  order  of  the  new  bond.  The  progress  of  the  concerted  reaction 

58  59 

is  then  described  as  ' the  sum  of  the  bond  orders  of  the  old 
and  the  new  bonds  is  maintained  (to  be  unity)  throughout  the 
rearrangement.  Such  an  interpretation  is  useful  but  not  essential 
and  alternative  views  of  n can  be  entertained. 

The  Interpolation  Formula 

The  thermodynamic  driving  force  for  the  reaction  is  the  free 
energy  difference  between  products  and  reagents.  The  simplest 
assumption  that  one  can  make  is  that  this  remains  true  also  during 
the  rearrangement  itself;  that  at  any  point  n enroute  from  reagents 
to  products  the  free  energy  G(n),  (measured  relative  to  the  reactants) 
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can  be  expressed  as 

G(n)  + G = (1-n) G + nG  . (7.1) 

r r p 

Here  the  subscripts  refer  to  reactants  and  products  respectively, 

G(n)  = 0 at  n=0  and  increases  linearly  with  n,  G(n)  = nAG,  reading 
AG  at  the  products,  n=l,  side. 

The  concept  of  a potential  energy  surface  provides  a theoretical 
molecular  interpretation  of  the  mechanism  for  the  concerted  loosen- 
ing of  the  old  bond  and  the  tightening  of  the  new  one/0,11  It 
also  shows  howtver  that  repulsion  between  the  non-bonded  atoms 
prevents  the  reaction  coordinate  from  running  entirely  down-hill 
in  the  exoergic  direction.  This  additional  source  of  inertia, 
equally  present  for  both  the  forward  and  reverse  reaction  requires 
that  the  linear  interpolation  (7.1)  be  augmented  by  a second  term, 
i.e. 

G(n)  = nAG°  + G?M(n) . (7.2) 

ci  • 

Here  the  n dependence  of  M(n)  is  to  be  such  that  G(n)  does  have  a 
maximum  at  some  value  of  n,  say  n=a.  The  magnitude  of  a (which, 
by  construction,  is  positive  and  between  zero  and  one)  is  the 
solution  of 

[3G(n)/anJn^a  =0,  (7.3) 

which,  using  (7.2),  is  seen  to  be  equivalent  to  the  previous 
definition  of  M(o) 

AG°/G°  = -l3M(n)/3n]n-a.  (7.4) 

As  AG°  varies,  so  will  the  magnitude  of  a which  is  determined  by 
(7.4).  To  insure  that  G(n)  is  always  maximal  at  n=ct  we  thus  need 

= G°32M/3n2  < 0. 

a 


to  require  that 
32G(n)/3n2 


(7.5) 
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The  condition  that  there  is  a maximal  free  energy  enroute  from 
reagents  to  products  is  this  sections  way  of  introducing  the 
inertia  of  the  chemical  reaction. 

The  Hammond  andLeffler  Postulates.  When  M(n)  is  independent  of 
AG°,  so  that  M(a)  depends  on  AG°  only  implicitly  (via  the  depend- 
ence of  a on  AG°)  it  follows  from  (7.2)  and  (7.4)  that 

[3G(n)/3AG°]n=a  -a.  (7.6) 

That  the  Br^nsted  coefficient  can  be  equated  to  the  magnitude  of 

the  progress  variable  at  the  point  of  maximal  free  energy  is 

12 

sometimes  called  theLeffler  postulate.  It  is  seen  to  be  a 

direct  consequence  of  the  form  (7.2)  for  G(n)  for  any  convex 

(cf.  (7.5))  function  M(n)  which  is  not  explicitly  a function  of 

AG°.  These  very  same  conditions  also  suffice  to  support  the  so- 

13 

called  'Hammond  postulate' . Specifically,  that  for  any  Br^nsted 
series,  the  location  of  the  point  of  maximal  free  energy  shifts 
towards  the  products  side  as  AG°  increases.  As  in  the  operational 
discussion,  (cf.  (6.15)),  the  proof  is  based  on  taking  the  deri- 
vative of  both  sides  of  (7.4)  with  respect  AG,  leading  to  3a/3AG° 

> 0.  The  advantage  of  the  model  is  that  it  identifies  the  Br^nsted 
slope  a,  which  satisfies  3a/3AG°>0,  with  the  location  of  the  point 
of  maximal  free  energy.  Any  irregularities  in  the  monotonic 
dependence  of  a on  AG°  are  thus  to  be  taken  as  evidence  that  one  is 


not  considering  reactions  of  the  same  series. 
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The  Potential  Energy  along  the  Reaction  Coordinate.  There  is  one 
important  limitation  to  the  discussion  based  on  (7.2) . The  opera- 
tional theory  requires  that  M(a)  depends  on  AG  only  via  its  a 
dependence.  This  is  insured  in  (7.2)  by  taking  M(n)  to  be  a 
function  only  of  n,  which  is  independent  of  AG.  This  assumption 
is  clearly  sufficient  to  prove  that  M(a)  is  a function  only  of  a 
but  it  is  not  necessary.  Of  course,  if  one  is  only  interested  in 
kinetic  behavior,  only  M(a)  is  needed,  and,  for  a given  series,  it 
is  only  a function  of  a.  If  however  one  wishes  to  regard  the 
static,  zero  temperature  limit,  then  G(n)  becomes  the  potential 
energy  as  a function  of  the  progress  variable  for  the  system 
(i.e.  the  potential  energy  along  a reaction  coordinate,  where  the 
coordinate  is  n) . It  then  becomes  a question  of  considerable 
interest  whether  it  is  possible  to  regard  M(n)  as  the  same  function 
for  the  entire  series  of  reactions.  An  answer  in  the  affirmative 
would  indicate  a very  useful  scaling  property  for  the  different 
potential  energy  surfaces  of  the  same  series.  There  are  at  present 
only  very  few  studies  related  to  this  point.  Starting  with  the 

BEBO  representation  for  the  energy  along  the  reaction  coordinate, 

17 

Marcus  has  shown  that  the  potential  energy  is  not  exactly  of 
the  form 

V (n)  = nA V+E°M ( n ) (7.7) 

Si 

with  M(n)  independent  of  the  difference,  AV,  in  bond  energies. 

However,  equation  (6.7)  is  satisfied  for  the  BEBO  form  of  M(n) 

and,  moreover,  the  AV  dependence  of  M(n)  is  at  best  small  and  can 

be  neglected^ ' ^ without  seriously  changing  V(n) . In  addition, 

28 

it  is  found  that  the  form  (7.7),  with  M(n)  being  only  n dependent 
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(as  in  (6.17))  can  represent  the  energy  profile  along  the  reaction 
coordinate  even  for  highly  exoergic  processes.  It  is  important  to 
note  in  this  respect  that  merely  examining  the  prediction  of  (7.7) 
for  the  correlation  of  the  barrier  height  (i.e.  V(a))  with  AV,  as  in 
figure  4,  will  not  suffice.  The  reason  being  that  the  result  that 
M(a)  is  only  a function  of  a defines  the  Br^nsted  series.  Figure  4 
is  not  a test  of  the  model  but  only  serves  to  determine  the  funct- 
ional dependence  of  M(a)  on  a.  Neither  is  the  proof  that  the 
Hammond  postulate  follows  upon  taking  the  derivative  of  M(a)  with 
respect  AG  and  the  same  applies  to  the  selectivity  reactivity  principle. 
All  these  results  are  valid  for  any  Br^nsted  series  of  reactions 
and  are  consistent  with,  but  do  not  establish  the  model. 


Energy  Consumption  and  Energy  Disposal.  The  model  enables  us  to 


correlate  the  location  of  the  free  energy  barrier  along  the  reaction 
coordinate  with  the  systematics  of  energy  requirements  and  energy 
disposal  in  a Br^nsted  series  of  reactions.  G =G  -G  has  been 
interpreted  in  section  II  as  the  excess  free  energy  of  the  reagents. 
Of  course,  G_  need  be  sufficient  to  overcome  AG°.  However,  even 

ci 

for  an  early  barrier  (a<l/2,  AG°<0) , G°  increases  monotonically  as 

cl 

a increases  and  reaches  G°  for  a symmetric  (a=l/2)  reaction.  The 

cl 

O 

excess  free  energy  in  the  nascent  products  is  G -AG  . It  is 
practically  negligible  for  a late  (a-*-l)  barrier  (highly  endothermic 
reactions  tend  to  produce  nearly  equilibrated  products) . As  the 
barrier  shifts  to  earlier  locations,  Ga~AG°  increases,  being  G° 
for  AG°=0  (at  a = 1/2) . 
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Figure  Legends 


1.  The  free  energy  of  activation  vs.  AG  , the  standard  free  energy 

change  in  the  reaction  for  three  different  Br^nsted  series  (i.e., 

for  three  different  values  of  G°  where  G°  = G at  AG°  = 0) . Since 

a a a 

the  Br^nsted  slope,  a = 3G  /8AG°,  is  independent  of  the  units  of 
energy,  G is  measured  in  the  same  units  as  AG°  and  the  units  in  the 

cL 

graph  need  not  be  specified  (changing  the  units  will  not  change 
the  graph  but  only  the  numerical  values  on  both  scales) . The  graphs 
are  drawn  for  the  particular  choice  M(a)  = -[alna  + (1  - a)ln(l  - a)]/ 
ln2  but  the  general  features  are  valid  for  any  convex  function  M(a) : 

Ga  increases  monotonically  with  AG°  with  a rate  of  increase  (the 
BrizJnsted  slope)  which  is  small  for  AG°  <<  0 and  increasing  as  AG° 
increases  (cf.  figure  2) . 

2 . The  dependence  of  the  Br^nsted  slope  on  AG°  for  the  three  series 
shown  in  figure  1.  As  is  evident  from  equation  (6.18),  a is  only 

a weak  function  of  AG°  for  AG°  values  in  the  range  |AG°|  < G°. 

Hence,  to  determine  a as  a function  of  AG  (and  hence  M(a))  from 
experimental  data,  it  is  necessary  to  measure  the  rates  both  for 
the  very  fast  (AG°  < -G°)  or  the  very  slow  (AG°  > G°)  members  of 

3 a 

the  series. 

3.  A Br^nsted  plot  (left  panel)  and  the  results  for  a vs.  AG° 

(right  panel)  for  the  X + H2  family  of  reactions  (X  a halogen  atom, 

T = 500K,  G°  = 4.08  kcal  mol  adapted  from  reference  28).  The 

cl 

52 

squares  are  the  measured  rates  for  the  X + H2  reactions  (a  slight 
extrapolation  of  the  Arrhenius  plot  for  F and  I was  required  to 
provide  results  at  a common  value  of  T for  the  entire  family) . The 
triangles  are  the  rates  for  XH  + H reactions  (obtained  using  detailed 
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balance)  . Strictly  speaking,  M(a)  should  be  determined  by  per- 
forming a Legendre  transform  of  the  G vs.  AG  plot  shown  in  the 

cl 

left  panel.  In  practice,  the  functional  form  (6.17)  for  M(a)  was 
assumed  and  only  the  value  of  G°  was  optimized,  to  yield  the  con- 
tinous  curve  shown.  The  Brjzfnsted  slope  (given  analytically  by 
(6.18))  is  shown  in  the  right  panel. 

4.  Variation  of  the  potential  energy  barrier  height  V(a)  with 
the  difference  AV  = AE°  in  bond  energies  for  the  X + Hj  series 
of  reactions.  (Energies  in  kcal  mol  1 , adapted  from  reference  25.) 
These  results  can  be  regarded  as  the  zero  temperature  limit  of 
figure  3.  The  continuous  curve  is  aAV  + V°M(a)  with  M(a)  and  a 

cl 

given  by  (6.17)  and  (6.18)  respectively  (V°  = 4.6  kcal  mol-1). 

a 

25 

The  points  shown  are  literature  values  for  V(a) . 
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Note  added  after  the  conference.  It  is  interesting  to  note  the 
connection  of  the  present  discussion  to  several  other  presentations. 
The  representation  of  the  exact  rate  constant  in  the  form  (1.2) 
has  also  been  discussed  by  Kuppermann.  The  concept  of  importance 
sampling  discussed  by  Muckerman  and  Faist  is  closely  related  to 
the  discussion  in  our  appendix  (although  the  purpose  of  the  two 
presentations  is  somewhat  different).  The  cycle  (3.11)  has  also 
been  considered  by  Weston  and  Klein.  Finally,  Dunning  has 
presented  a correlation  similar  to  Figure  4 except  that  his  barrier 
heights  are  the  results  of  accurate  ab  initio  computations  whereas 
the  points  shown  in  Figure  4 (and  in  reference  25)  are  semi-empirical 
estimates  which  best  account  for  the  experimental  results. 
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Appendix 

The  appendix  provides  a derivation  of  the  result  kr (T ) = 
kTST(T)K(T)  where  < (T)  , (k  (T)  <1)  is  the  transmission  coefficients 
along  the  same  lines  used  to  derive  the  result  kr(T)=Zpr(T)  in 
section  V.  As  will  become  clear,  the  two  results  indicated  are 
only  two  special  cases  of  an  infinite  number  of  expressions  all  of 
which  share  the  common  form,  that  of  an  upper  bound  to  the  reaction 
rate  constant  times  a reactivity  (or  a transmission)  factor.  There 
are,  of  course,  particular  reasons  for  centering  attention  on  the 
two  special  forms  (Z  is  independent  of  AG;  < (t)  is  chosen  to  be  as 
nearly  unity  as  possible)  but,  in  principle,  other  choices  are 
possible. 

The  derivation  starts  with  the  general  form  (cf.  (5.8)) 

2°° 

k(n,T)=irB  /dE.  f (E  ) v[N  (n,E  )/N(n,E.)]  (A.l) 

0 t t r t ^ 

With  the  sampling  distribution 

N(n,Efc)  = f(Efc)N(n)  (A. 2) 

one  obtains 

k(n,T)  = 7rB2Jr  (n)/N(n)  (A. 3) 

where  Jr (n) , 

00 

J_ (n)  = / vN  (n,E  ) dE  (A. 4) 

r 0 r 

is  the  flux  of  reactive  reactants  when  the  reactants  are  sampled 
as  in  (A. 2).  If  Jfa(n)  is  an  upper  bound  to  Jr(n), 

kr(n,T)  = kfa(n,t) [Jr(n)/Jb(n)l  (A. 5) 

Here  kb(n,T3is  defined  as  in  (A. 3)  with  Jb(n)  replacing  Jr (n) 
and  the  term  in  the  square  brackets  is,  by  construction,  bounded 


I 

r 
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by  unity.  Proceeding  to  the  thermal  rate  constant 

kr(T)  = £[N(n)/N]kr(n,T) 
n 

= TrB2Jr(T)/N 
= kb(T)  [Jr(T)/Jb(T)J 


(A. 6) 


where 


and 


Jr(T)  = Z[N(n)/N]Jr(n,T) 
n 


Jb(T'  = 1 [N(n)/N] Jfa(n,T)  > Jr(T) 
n 


(A. 7) 


(A. 8) 


(A. 9) 


Since,  by  definition 

kb(T)  = irB2Jb(T)/N 
it  follows  from  (A. 8)  and  (A. 6)  that 

kb(T)  = E[N(n)/N]kb(n,T)  . (A.  10) 

Two  particular  choices  for  Jb(T)  are  (i)  the  flux  of  all 
trajectories  that  originate  from  the  reactants.  Then  kb(T)  = Z, 
leading  to  the  results  of  section  V.  (ii)  the  flux  of  all  traject- 
ories that  originate  from  the  reactants  and  that  proceed  to  cross 

. 47-49 

a given  surface  on  route  from  reagents  to  products  . This 

defines  J+(T)  such  that  k+(T)  is  the  transition  state  rate  constant 
and  the  transmission  factor  is  Jr (T) /J+ (T) . If  there  is  a surface  of 
no  return  than  for  that  surface  Jr(T)  = J+(T.)  and  <(T)  = 1.  Other- 
wise, k(T)«  1.  The  free  energy  of  activation,  G+(T)-G(T),  of  the 
transition  state  theory  rate  constant  is  thus  a lower  bound  to  the 
correct  value  G (T)  , and  equals  G (T)  when  ic(T)  = 1.  The  criterion 

a 6L 

of  locating  the  surface  of  no  return  by  maximising  G+(t)  has  indeed 
been  extensively  discussed*®' 

Finally,  it  is  important  to  note  that  for  a bimolecular  reaction 


49 


the  free  energy  of  activation  G is  independent  of  the  choice  of 

concentration  units  for  the  rate  constant.  In  other  words,  Gr 

(where  G =G  -G)  is  closely  related  to  but  not  identical  to  G^ 
cL  r 

defined  by3'9'15 

kr(T)  = (kT/h) exp (-G^/RT) /Q  (A. 11) 

where  Q is  the  partition  function  for  the  reactants.  To  establish 
the  relation  it  is  convenient  to  use  the  identity9 

(kT/h)QrR  = ZQrt  (A. 12) 

Here  Q . and  Q D are  the  translational  and  rotational  partition 

ru  iK 

functions  for  the  reactive  reactants  (where  the  reduced  mass  for 
relative  motion  is  identical  for  the  reactants  and  reactive 
reactants).  Comparing  (A. 11)  with  kr(T)  = Zexp  (-Gr/RT) /Q  and  using 
(A. 12) 

Gr  + RTlnQrt  H Gr-Grt  (A. 13) 

= G+  + RTlnQrR  = G+-GrR 
or 

- Gr-Grt+GrR  (A. 14) 

For  the  reactive  reactants  the  motion  is  separable  into  translational 
and  internal  contributions.  Hence  one  can  factor  Qr  as  QrtQrI  where 
QrI  is  the  internal  partition  function.  It  follows  that  Gr-Grt 
(cf.  (A. 13))  is  the  internal  free  energy  of  the  reactive  reactants. 
For  the  transition  state  the  separable  approximation  is  not  necessar- 
ily valid49' 67 * If  it  is,  one  can  interpret  G^-G  D as  the 

JTK 

internal  free  energy  at  the  transition  state.  One  should  note 

35 

however  that  one  can  always  express  the  rate  constant  as  in  (A. 11), 
irrespective  of  the  validity  of  transition  state  theory,  so  that 


(A. 14)  can  be  taken  as  a definition  of  G^.  The  advantage  of 

G —G  -G  over  AG^=G^-G  is  that  the  standard  state  dependence 
a r 

of  G t is  exactly  cancelled  by  a similar  dependence  of  the  Gt 
contribution  to  G,  so  that  G„  is  independent  of  concentration  units. 
Not  so  for  AG^.  The  elimination  of  Grt  (cf.  (A. 14))  means  that 

+ CO  CQ 

has  the  dimensions  of  ' energy.  [ dimensionless  factor  + log 


(volume) ] . 
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Abstract 

The  vibrational  excitation  of  an  anharmonic  (Morse)  oscillator  in  a collinear 

collision  is  determined  analytically  using  a dynamical  algebra  consisting  of  three 

operators.  Only  one  parameter  is  however  required  to  specify  the  transition 

probabilities.  The  dependence  of  V-*T  processes  on  the  initial  vibrational  state 

is  examined.  A weak  coupling  scaling  law  P , , * n[l-(n-l) /2k]P, where  k is  the 

rrm-i  l-H) 

number  of  bound  states  is  derived.  At  finite  coupling  strengths  the  decline  of 
Pn-+n-l^n  wit^  n is  faster  than  that  suggested  by  the  scaling  law. 
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I.  Introduction 


Transition  probabilities  in  molecular  collisions  often  show  considerable 
deviance  from  that  expected  if  the  final  phase  space  is  uniformly  populated. 
Surprisal  analysis  [1]  has  demonstrated  that  such  deviances  can  often  be  charact- 
erized as  due  to  simple  constraints.  Recently,  it  was  shown  that  such  constraints 
can  be  identified  on  the  basis  of  exact  dynamical  considerations  [2,3].  A minimal 
set  of  constraints,  the  'dynamical  algebra',  was  found  sufficient  to  determine  the 
scattering  matrix  [3],  Brute  force  attempts  to  construct  the  dynamical  algebra 
for  realistic  potentials  often  become  cumbersome.  The  purpose  of  this  letter  is 
to  draw  attention  to  a powerful  systematic  method  for  the  construction  of  such 
algebras  [4-6]  and  to  demonstrate  its  application  up  to  and  including  drawing 
conclusions  of  physical  interest. 

Consider  a collinear  collision  of  a particle  with  a Morse  oscillator.  In 
the  classical  path  approximation  this  collision  can  be  described  as  the  evolution 
of  a Morse  oscillator  under  a time  dependent  force  (which  vanishes  at  t -*•  ±°°) . We 
obtain  an  analytic  expression  for  the  transition  probabilities  and  compare  with  the 
corresponding  results  for  a harmonic  oscillator.  The  comparison  is  carried  out  in 
terms  of  the  'reduced  anharmoniclty'  or  maximal  number  of  bound  states 

k - (2  pDo)1/2/ah 

which  arises  as  a natural  parameter  in  the  formulation  of  the  problem.  (Here  Dq 
is  the  ground  state  dissociation  energy,  p the  oscillator  mass  and  a is  the  range 
of  the  Morse  potential  [7]).  As  k -*•  00 , the  results  reduce  to  those  [3,8]  of  the 
forced  harmonic  oscillator.  The  most  significant  effect  for  a finite  k is  the 
modification  of  the  familiar  weak-coupling  scaling  law  of  V ♦ T transfer.  Here 

Vn-l  ' ■[W-1),!l|,W 

where  n is  the  initial  vibrational  state. 

In  addition  to  obtaining  explicit  formulas  for  the  quantal  probabilities,  we 
also  obtain  their  classical  limit  and  demonstrate  that  there  is  a limited  band 
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of  final  states  which  are  allowed  by  classical  mechanics.  The  width  of  the  band 
increases  with  increasing  coupling  strength,  up  to  the  full  width  allowed  by  quantal 
dynamics . 

II.  Scattering  Theory 

The  explicit  construction  of  the  scattering  matrix  in  terms  of  the  operators  in 
the  dynamical  algebra  has  been  presented  elsewhere  [3].  Here  we  only  state  the  final 
result  and  comment  on  its  adaptation  to  those  dynamical  algebras  determined  by  the 
systematic  procedure  [4-6],  A detailed  presentation  of  the  systematic  procedure  for 
the  generation  of  dynamical  algebras  with  special  reference  to  scattering-theoretic 
applications  is  in  preparation. 

For  a collision  problem  in  the  classical  path  approximation,  consider  the 
following  form  for  the  Hamiltonian 

H = E h (t)H  (1) 

_ n n 
n»0 

Here  the  Hn's  are  operators,  with  H^  being  the  Hamiltonian  of  the  unperturbed 

system.  The  hn(t)'s  are  function  of  time,  which,  for  ni*0,  vanish  both  before  and 

after  the  collision.  A dynamical  algebra  is  a set  of  operators  H which  is  closed 

n 

under  commutation,  [3] 

[H  ,H  ] - E Cr  H , (2) 

1 n’  mJ  r nm  r* 

and  is  therefore  a Lie  algebra  [9].  If  necessary,  the  set  of  operators  that  appear 
in  H can  be  augmented  by  the  introduction  of  additional  operators,  until  the  closure 
condition  (2)  obtains.  For  these  additional  operators,  the  h(t)  coefficients  in 
(1)  are  to  be  taken  as  zero  throughout. 

Given  a dynamical  algebra,  the  evolution  operator  is  of  the  form  (equation  $.3) 
of  ref.  3) 

U(t,  t') -exp  [-(i/ll)  Eun(t,t')Hn]  (3) 

and  the  scattering  matrix  is  obtained,  as  usual,  by  evaluating  the  matrix  elements 
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of  U (in  the  interaction  picture)  and  taking  the  limits  t-*»,  t'-*-00).  Explicit 
equations  of  motion  for  the  coefficients  u^t.t')  have  been  derived  (equation  (5.19) 
of  ref.  3). 

A systematic  procedure  for  the  construction  of  a Lie  algebra  for  a given 

Hamiltonian  has  been  presented  [4-6],  In  particular,  it  is  possible  to  transform  the 

original  Schrodinger  operator  K=H-i9/3t  to  R 

R = DKD-1  (4) 

such  that  the  algebra  for  K is  far  simpler.  This  is  particularly  so  if  the  spectrum 

of  Hq  is  known.  (See  the  examples  of  ref.  5).  One  can  then  determine  the  evolution 

operator  for  R and  transform  back  to  obtain  the  evolution  operator  for  K.  In 

practice  one  is  only  interested  in  matrix  elements.  Specifically,  if  ip  (t)  is  a 

n 

solution  of  the  time  dependent  Schrodinger  equation  for  Hq,  we  need  the  matrix 

elements  of  U(t,t')  between  the  states  i(/  (t)  and  ip  (t').  But  these  are  the  same 

n m 

as  the  matrix  elements  of  0(t,t')  in  the  transformed  basis  D4>n(t).  (We  require 
-1 

that  DHgD  = Hq) . Solving  the  scattering  problem  for  K provides  therefore  the 
solution  of  the  scattering  problem  for  K.  All  that  is  required  for  the  original 
problem  is  to  determine  the  S matrix  in  the  transformed  basis. 


III.  The  Forced  Horse  Oscillator 

The  bound  state  energy  levels  of  the  Morse  oscillator  are  17] 

En  - -(aV/2y)(k-n)2  . (5) 

k is  a convenient  reduced  measure  of  the  anharmonlcity  and  we  shall  label  the 
bound  states  by  both  k and  n,  although  k is  fixed  for  any  given  physical  applicat- 
ion. The  solution  of  the  time  dependent  Schrodinger  equation  for  an  unperturbed 

oscillator  will  thus  be  written  as  ip  (t)  = exp(-iE  t^H)|n,k>. 

n n 

The  similarity  transformation  employed  here  is  that  of  time  dilation  [5,6] 
such  that 

D^n(t)  - exp[-i(n-k)t] |n,k>.  (6) 
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The  transformation  linearises  the  spectrum  of  H^. 

The  forcing ' of  the  oscillator  is  described  by  operators  which  change  its 
state.  In  general,  the  oscillator  can  also  dissociate.  To  simplify  the  discussion 
we  shall  use  an  algebra  that  excludes  this  possibility  and  introduce  an  auxiliary 


label  m such  that 

n = k-|m|  (7) 

so  that  m = ±k  is  the  ground  state  and  m takes  2k+l  values  from  -k  to  +k.  The 

forcing  terms  will  induce  Am  transitions, 

Q+|k,m>=c+(km) |k,m±l>  (8) 

Here  c+(km)  is  a numerical  coefficient.  To  insure  that  the  problem  reduces  to  that 

4* 

of  the  forced  harmonic  oscillator  and  that  Q+  = Q_  we  take 

c±(km)  = [k(k+l)-m(m±l)]1/2  (9) 

It  follows  from  (8)  and  (9)  that 

[Q+»Q_]|k,m  > = 2m | k , m>  (10) 

Hence  our  algebra  also  contains  the  operator  Qq, 

Qo|k,m>=m|k,m>  (11) 

or 

[Q+.Q-]  - 2Qq.  (12) 

Are  there  additional  operators?  By  examining  the  action  of  [Qq,Q+]  on  the  basis 
functions  one  readily  verifies  that 

[Q0,Q±]  = ±Q±  (13) 


An  algebra  for  the  transformed  Hamiltonian  is  given  by  {Q+*Q_»Qq}  and  is  that  of 
SO (3).  An  explicit  coordinate  space  realization  shows  that  Q+  are  linear  combinations 
of  the  oscillator  coordinate  and  momenta  with  coefficients  that  depend  on  the 


coordinate.  In  the  limit  k ■*  00  they  reduce  to  the  algebra  of  the  forced  harmonic 
oscillator  [3].  In  that  limit 
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where  a is  the  harmonic  oscillator  raising  operator  and  similarly  for  Q . Hence 
[Q+.QJ  ->-2k[a;a]  = 2kl  (15) 

so  that  Qq  -►  kl  where  I is  the  identity  operator.  Now 
<m+l,k| [Qq,Q+] |m,k>  = c+(k,m) 

21/2(k)3/2<n+l| [I,a+] |n>  = k1^2  ( 6) 

or  [I,a^]  = 0,  QED.  The  algebra  we  have  generated  does  reduce,  as  k 00  to  the 
dynamical  algebra  {a^,a,l}  of  the  forced  harmonic  oscillator. 


IV.  Transition  Probabilities 

The  determination  of  the  transition  probabilities  requires  the  representation 
of  the  evolution  operator  in  the  basis  Di/>n(t).  The  most  general  form  of  the  scatter- 
ing operator  for  the  dynamical  algebra  (Q+,Q_,Q0)  is 

S « exp[i(qQ+  + q*Q_  + 6Qq) ] (17) 

where,  since  (Q+)^  = Q_,  S is  unitary.  One  can  rewrite  S in  a product  form  (cf. 
equations  (4.4)  or  (5.3)  of  [3]),  and  its  most  general  form  is  ([10],  pages  54-55 
equation  (4.1.8)  in  particular) 

S = exp[iaQo]exp[8(Q+-Q_)]exp[iyQo]  (18) 

whose  representation  is  [10] 

<m'k|s|mk>  = (aBy)  = exp(im'*y)d^  (8)exp(ima)  (19) 

mm  mm 

Here  the  D's  are  the  unitary  rotation  matrices  and  the  d's  are  the  reduced  unitary 

rotation  matrices.  The  symmetry  properties 

d (6)  - (-l)m""md^(B)  = (-l)m'"md(H  (6)  (20) 

m m mm  -m  — m 

will  be  invoked  below. 

The  transition  probabilities  among  different  oscillator  states  are  given  by  the 
squares  of  the  S matrix  elements 

T>  — I A (k)  _ 


k-„-  <»i 


,-2k,_  (-l,V2v~n'H>  1/2 


nt  (2k-ft)  ! (o')  ! (2k-n')  t v7(v-n4ti't) ! (n-v)  I C2k-n',-v)  I 
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where  p = tn  (8/2) . They  are  symmetric  (cf.(20))  and  normalised. 

There  is  only  one  group  parameter  (i.e.  $)  that  enters  the  expression  for  the 

transition  probabilities.  In  principle,  it  can  be  computed  from  the  equations  of 

motion  for  the  group  parameters,  [3],  in  terms  of  the  forcing  functions  h+(t)  (cf. 

equation  (1)).  For  the  purpose  of  interpreting  the  results  it  is  sufficient  however 

to  note  that  to  lowest  order  in  the  forcing  terms  (cf.  (18)) 

SDW  = exp[iaQo][I+B(Q+-Q_)]exp[iYQo]  (22) 

It  follows  that  in  first  order  the  perturbation  induces  transitions  only  among  adjacent 

2 

states  of  the  oscillator  and  that  (cf.  (16))  the  proper  coupling  parameter  is  k8  . 

The  final  vibrational  state  distribution  is  monotonic  for  transitions  out  of  the 

ground  state  (cf.  (25)  below)  but  shows  a typical  quantum  interference  for  n / 0,  as 

is  shown  in  Figure  1.  In  the  classical  limit,  these  oscillations  can  be  averaged 

over  [11].  The  vector  model.  Figure  2,  gives  the  classical  distribution  as 

P ^ = (irk)  ■*"[  (cosS-cosvcosvO^-sin^vsin^v']  (23) 

n-*n 

where  cosv=l-n/k.  The  most  probable  final  states  in  the  classical  limit  are  the  two 

roots  of  8=v'±v,  which  define  the  end  points  in  the  vector  model,  figure  2. 

In  the  k (or  harmonic  oscillator)  limit,  we  obtain  (cf.  (14)) 

S = exp[a(a+-a)]  (24) 

1/2 

witha=(2k)  8,  which  is  the  previously  derived  result  ([3],  equation  (3.40)). 
Explicitly,  and  using  the  0 -*■  n transitions,  we  have  obtained  here  (cf.  (21)) 

P0  ^ n=  (2n)[pn/(1+P)2k]*  (25) 

Using  the  Stirling  approximation  this  can  be  shown  to  reduce  to 

PQ  >a  = [gn/n!  ]exp(-g) , k -*•  ® (26) 

where  g = 2kp(-*-2k$  as  fW)) , which  is  the  well  known  result  [8]  for  the  forced 

harmonic  oscillator. 

Explicit  expressions  for  the  mean  final  quantum  number  (or  higher  moments 
thereof)  are  readily  obtained  by  differentiating  the  normalisation  condition 


EnP0+n 


1 with  respect  to  p.  One  then  finds 
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V.  Discussion 

The  role  of  the  anharmonicity  is  best  demonstrated  by  comparing  the  results  for 

the  harmonic  and  the  Morse  oscillators.  Thus  the  ratio  of  the  0*-l  transition 

2k 

probabilities  (given,  using  (25)  and  (26)  by  (1+p)  exp(-g)),  is  always  below  unity. 


A particularly  interesting  ratio  is  that  of  the  de-excitation  probabilities 
P | q * For  the  collinearly  forced  harmonic  oscillator  and  to  the  lowest  order 

in  the  coupling  strength  this  ratio  is  [12]n.  For  this  reason  it  has  become  customary 
[13]  to  scale  the  n-+n-l  transition  probabilities  by  n.  It  is  thus  of  interest  to 

I 

note  that  the  same  ratio  for  the  Morse  oscillator  (in  the  weak  coupling  limit)  is 
given  by 

| c_ (km) | 2 f | c_ (kk) | 2-n [ 1- (n-1) /2k]  (29) 

or  P_ . ,/n  is  a decreasing  function  of  n.  Of  course  the  modified  anharmonic 

rrm-i 

scaling,  ^ | ^_^/n[l-(n-l) /2k]  is  only  valid  in  the  weak  coupling  limit.  Even  so 

it  does  account  for  the  often  observed  decline  of  P ,/n  with  n,  and  does  compare 

n-+n-l 

favorably  with  the  trends  noted.  Figure  3,  using  transition  probabilities  computed 
using  the  full  form  (21) . The  corresponding  result  for  excitation  is 

(30> 

•k 

An  improved  scaling  which  remains  valid  not  only  as  p-*0  but  as  long  as  nkp<l  is 

*In  applying  these  scaling  laws  it  should  however  be  noted  that  8 (cf.  equation 
22)  was  taken  to  be  the  same  for  the  1«-K)  and  the  n<  >n±l  transitions.  This  is 
correct  in  the  classical  path  approximation.  However,  when  the  relative  motion  is 
treated  quantum  mechanically  one  expects  8 to  increase  somewhat  as  n increases 
(reflecting  the  smaller  energy  gap  for  the  higher  transitions.  Approximate  procedures 
for  incorporating  this  aspect  in  the  classical  path  formulation  have  been  proposed  [141). 
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Pn-^-l“nPl-+Ol  (2k'n"1)/2k1 1 1_(2k"n)  (n'1)  <P/2>  l2  (31) 

obtained  by  retaining  the  first  two  terns  in  the  sum  in  (21) . 

VI . Summary 

In  demonstrating  the  dynamical  significance  of  the  information  theoretical 
approach  it  was  argued  [3]  that  one  can  determine  the  functional  form  of  the 
scattering  matrix  (or  of  the  distribution  of  final  states)  without  solving  the 
equations  of  motion.  The  resulting  scattering  matrix  is  then  obtained  as  a function 
of  the  group  parameters  (here  a, 6 and  y»  cf.  equation  (18)).  The  solution  of  the 
equations  of  motion  is  required  only  for  the  determination  of  the  numerical  values 
of  the  group  parameters  [3]. 

In  this  paper  an  explicit  construction  of  the  scattering  matrix  and  an  examin- 
ation of  the  resulting  transition  probabilities  was  carried  out  for  a non-trivial 
collision  problem.  A key  ingredient  has  been  the  application  of  a systematic 
procedure  for  the  generation  of  Lie  algebras  [4-6].  In  particular,  one  can  construct 
dynamical  algebras  that  do  allow  the  Morse  oscillator  to  dissociate  or  for  other 
types  of  anharmonic  oscillators.  A complete  account  of  the  scattering  theoretic 
implications  of  the  systematic  procedure  for  the  construction  of  Lie  algebras 
associated  with  a given  Hamiltonian  is  in  preparation. 
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Figure  Legends 


1.  Transition  probabilities  for  the  forced  Morse  oscillator  (at  the  realistic 
value  k * 15  and  intermediate  coupling  strength  6 * 36.9°).  The  destructive 
quantal  interference  for  n * 3 is  evident.  The  classical  probabilities 
(equation  (23))  are  shown  for  comparison.  The  states  n = 0 and  n > 5 are 
classically  forbidden. 

2.  The  vector  model  for  the  classical  distribution  of  final  vibrational  states, 

based  on  the  classical  limit  of  the  d matrices  [11].  ON  and  ON'  are  both  of 

length  k.  The  classical  P ..  distribution  is  obtained  by  projecting  the  tip 

irn 

of  the  k vector  onto  the  ON'  axis,  while  the  vector  k spans  a circle  (perpen- 
dicular to  the  NON'  plane)  centered  at  a distance  n from  N.  The  classical 
turning  points,  where  the  probability  is  maximal,  are  at  v'±v*{5,  and  the  range 
of  n'  values  allowed  in  the  classical  limit  are  indicated.  The  increase  of  the 
range  of  higher  vibrational  states  that  are  allowed  by  classical  mechanics 
clearly  increases  as  the  strength  of  the  coupling  (<*&)  increases.  The  transition 
to  the  highest  bound  vibrational  state  (n'=k)  is  classically  allowed  only  for 
B^90-v . 

3.  The  scaling  law  for  n+n-1  transitions  for  k * 15.  In  the  harmonic  oscillator 
weak  coupling  limit  all  the  points  should  be  at  unity.  Solid  line:  the  Morse 
oscillator  weak  coupling  limit.  For  finite  but  weak  (p  ■ 1/500  or  g * 0.06 
and  p * 1/150  or  g = 0.1)  couplings,  the  decline  is  still  monotonic  and  is  well 
represented  by  (31).  At  intermediate  coupling  (p  - 1/60  or  g - 0.5),  the 
behavior  is  no  longer  monotonic. 
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